Mon. Not. R. Asteon. Soc. 0Q0.mi38l(2008) Printed 27 January 2009 (MN MgC style file v2.2) 



Stellar Population and Kinematic Profiles In Spiral Bulges & Disks: 
Population Synthesis of Integrated Spectra 

Lauren A. MacArthur 1 *, J. Jesus Gonzalez 2 f, and Stephane Courteau 3 ^ 

^Department of Astrophysics, California Institute of Technology, MS 105-24, Pasadena, CA 91125 

2 Instituto de Astronomia, Universidad Nacional Autonoma de Mexico, Apdo Postal 70-264, Cd. Universitaria, 04510 Mexico 
3 Department of Physics, Engineering Physics & Astronomy, Queen 's University, Kingston, ON K7L 3N6, Canada 



Accepted Jan 16, 2009 



ABSTRACT 

We present a detailed study of the stellar populations (SPs) and kinematics of the bulge 
and inner disk regions of eight nearby spiral galaxies (Sa-Sd) based on deep Gemini/GMOS 
data. The long-slit spectra extend to 1-2 disk scale lengths with S/N/A ^ 50. Several differ- 
ent model fitting techniques involving absorption-line indices and full spectrum fitting are 
explored and found to weigh age, metallicity, and abundance ratios differently. We find that 
SPs of spiral galaxies are not well matched by single episodes of star formation; more rep- 
resentative SPs must involve average SP values integrated over the star formation history 
(SFH) of the galaxy. Our "full population synthesis" method is an optimised linear combi- 
nation of model templates to the full spectrum with masking of regions poorly represented 
by the models. Realistic determinations of the SP parameters and kinematics (rotation and 
velocity dispersion) also rely on careful attention to data/model matching (resolution and flux 
calibration). The population fits reveal a wide range of age and metallicity gradients (from 
negative to positive) in the bulge, allowing for diverse formation mechanisms. The observed 
positive age gradients within the effective radius of some late-type bulges helps reconcile the 
long-standing conundrum of the co-existence of secular-like kinematics, light profile shape, 
and stellar bar with the "classical"-like old and a-enhanced SPs in the Milky Way bulge. 
The disks, on the other hand, almost always show mildly decreasing to flat profiles in both 
age and metallicity, consistent with inside-out formation. Our spiral bulges follow the same 
correlations of increasing light-weighted age and metallicity with central velocity dispersion 
as those of elliptical galaxies and early-type bulges found in other studies, but when SFHs 
more complex and realistic than a single burst are invoked, the trend with age is shallower and 
the scatter much reduced. In a mass-weighted context, however, all bulges are predominantly 
composed of old and metal-rich SPs. While secular contributions to the evolution of many of 
our bulges is clearly evident, with young (0.001-1 Gyr) SPs contributing as much as 90% of 
the optical (V-band) light, the bulge mass fraction from young stars is small (< 25%). This 
implies a bulge formation dominated by early processes that are common to all spheroids, 
whether they currently reside in disks or not. While monolithic collapse cannot be ruled out 
in some cases, merging must be invoked to explain the SP gradients in most bulges. Further 
bulge growth via secular processes or "rejuvenated" star formation generally contributes min- 
imally to the stellar mass budget, with the relative secular weight increasing with decreasing 
central velocity dispersion. 

Key words: galaxies: spiral — galaxies: evolution — galaxies: formation — galaxies: stellar 
content — galaxies: bulges 



1 INTRODUCTION 
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In the context of the currently favored cosmological A-Cold Dark 
Matter model (ACDM) of our Universe, the formation and evolu- 
tion of galaxies remains a major unsolved problem. In this model, 
structures form in a bottom-up, hierarchical, manner by which 
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smaller fragments merge together to form more massive systems. 
While this model has been very successful at reproducing observa- 
tions on large scales, a number of issues remain to be resolved at 
the galaxy-scale regime (Moore et al. 1999; Navarro & Steinmetz 
2000; Primack 2007, and references therein). In particular, the for- 
mation of disk galaxies is not well represented in current imple- 
mentations of simulations based on ACDM (e.g. Bell et al. 2003; 
Kaufmann et al. 2007; Dutton et al. 2007). The discrepancies do 
not necessarily indicate a failure in the ACDM model, but rather 
point to physical regimes and processes that are either not well un- 
derstood or difficult to implement in large simulations. Given the 
theoretical difficulties and limitations faced by current models, a 
true understanding of disk galaxies must include information from 
an observational perspective. 

All aspects of baryonic physics, from hydrodynamical pro- 
cesses involving gas, through star formation and feedback from su- 
pernovae and AGN, can play a role in the formation and regulation 
processes that shape galaxies, particularly where dynamically vi- 
olent processes such as major mergers and rapid collapse are not 
dominant. For spiral galaxies, the very presence of a disk compo- 
nent implies that no such major event has occurred since the forma- 
tion of the dynamically fragile disk. Additionally, most disk galax- 
ies harbor central bulge components whose prominence in terms 
of total mass or light weight spans a wide range from truly bulge- 
less systems to the bulge-dominated early-type SO/Sa's. A question 
that naturally arises is if the final step to the disk-less, i.e. pure el- 
liptical, galaxies is a natural extension of bulge-to-total ratios in 
systems along the Hubble sequence. Such an extension would im- 
ply that the dominant bulges of early-type spirals have followed a 
similar evolutionary path to those of pure elliptical galaxies. How- 
ever, coming from the opposite end, a connection between the tiny 
bulges of late-type spirals, through early-types, to pure ellipticals 
may not seem as obvious, and allowance must still be made for 
the pure disk galaxies in any scenario. The emerging observational 
picture of bulge formation currently implies similar evolutionary 
paths for early-type spirals and pure ellipticals, while the later-type 
bulges exhibit much more diversity in their observed properties, 
consistent with having formed through a secular redistribution of 
material from the disk (e.g. Kormendy & Kennicutt 2004). The 
nomenclature predominant in the literature to distinguish between 
formation scenarios for bulges names those appearing very similar 
to pure elliptical systems (i.e. rapid and/or violent formation which 
includes both the monolithic collapse and major merging scenarios) 
as "classical bulges" and those formed secularly from the disk as 
"pseudobulges". However, uncertainties in making a clear-cut dis- 
tinction between these cases still remain, particularly in regards to 
the stellar populations (SPs) of spiral bulges. 

A very useful probe in discerning between formation scenarios 
is a detailed breakdown of the age, metallicity (Z), and kinematic 
properties of the SPs comprising bulges of all types. Information 
about both light- and mass-weighted quantities is needed to form 
a comprehensive picture of the star formation history (SFH) of a 
given system. However, beyond our Local Group, we are limited 
to observations of the integrated light along a given line of sight. 
Of relevance is whether the integrated light at any location in a 
galaxy can be deconvolved into the relative fractions of stars of 
a given population that contribute to the total light. This question 
is especially acute for spiral galaxies that are known to harbor a 
mixture of young and old stars and can suffer from the extinction 
effects of dust. 

Observations of colour gradients and absorption line-indices 
in galaxies can both be interpreted as a single-burst, single- Z stel- 



lar population (SSP), from which light-weighted SSP-equivalent 
ages and metallicities can be drawn with the help of spectral syn- 
thesis models (e.g. Vazdekis 1999; Bruzual & Chariot 2003, here- 
after BC03; Thomas, Maraston, & Bender 2003; Le Borgne et al. 
2004; Maraston 2005; Schiavon 2007). Elliptical galaxies (e.g. 
Trager, Faber, & Dressier 2008) and globular clusters (e.g. Cohen, 
Blakeslee & Rhyzhov 1998; Puzia et al. 2005; Hernandez & Valls- 
Gabaud 2008), are often assumed to result from the evolution of a 
single burst of star formation and thus can be represented by a SSP 
of a given age and Z. This assumption, however, breaks down for 
spiral galaxies that are believed to have been converting gas into 
stars at a relatively constant rate over most of their lifetime (e.g. 
Kennicutt 1983; James, Prescott, & Baldry 2008). Since the pres- 
ence of a young SP will dominate the optical light even when its 
contribution to the mass is minimal, the light-weighted ages will 
not be representative of the entire population. Even pure ellipticals 
show signs of "frostings" of young SPs and may not be well repre- 
sented by a single SSP (Serra & Trager 2007). 

Greater potential for disentangling the SFH of spiral galax- 
ies exists if they can be assumed to have proceeded smoothly with 
time. This approach was followed by MacArthur et al. (2004) us- 
ing broad-band optical colours for 172 low-inclination disk galax- 
ies. To the extent that the nearby spiral galaxies in this sample 
share a similar underlying SFH and that no burst involving more 
than ~ 10% of the galaxy mass has occurred within the past 1- 
2Gyr, the method provides reliable relative results. Among others, 
it was found that the SFH of spiral galaxies depends strongly on the 
galaxy potential and halo spin parameter. 

Conclusions based on colour gradients are, however, plagued 
by a degeneracy between the effects of age, Z, and dust, all leading 
to redder colours. Spectroscopic techniques, such as those based 
on absorption-line equivalent widths that are largely impervious to 
dust effects (MacArthur 2005), or fitting of the full spectral energy 
distribution (SED) which can incorporate a dust component, offer 
a more detailed and discriminating view, especially in light of the 
latest implementations of SP synthesis models (e.g. Vazdekis 1999; 
BC03; Le Borgne et al. 2004). 

Early-type spirals with their large bulge-to-disk ratio offer 
a significant observational advantage over late-types since their 
bulges, which rise significantly above and below the disk, can be 
studied spectroscopically in the edge-on perspective free from disk 
contamination and extinction from dust. Late-type bulges, on the 
other hand, are small and often no thicker than the galaxy disk it- 
self, thus requiring a face-on projection for their study. While dust 
and contamination from the inner disk will still thwart any pristine 
observation of the bulge, the face-on orientation minimizes line-of- 
sight integrations and keeps the light distribution of the bulge and 
disk free of inclination effects. 

Spectroscopic studies of bulges spanning the full range of 
Hubble types are few, and the results are often conflicting. The 
absorption-line studies of the central regions of Sa-Sc spirals of 
Trager et al. (1999) and Proctor & Sansom (2002, hereafter PS02) 
both find that late-type bulges cannot be reproduced using primor- 
dial collapse models and invoke extended gas infall onto the cen- 
tral bulge to explain the observations. In a similar study, Goudfrooij 
et al. (1999) conclude the opposite, their findings being more com- 
patible with predictions of the 'dissipative collapse' model than 
with those of the 'secular evolution' model for bulge formation. 
Their sample is, however, dominated by early-type spirals and these 
conclusions may not apply to the full range of spiral types. Finally, 
the analysis of Moorthy & Holtzman (2006) of line strengths in the 
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bulges and inner disks of 38 spirals also favours a scenario whereby 
merging is the dominant mechanism for bulge formation. 

Thomas & Davies (2006) reanalyzed the PS02 data and found 
no difference between the SPs of spiral bulges and Es at a given 
central velocity dispersion, ao, concluding that processes involving 
disk material cannot be responsible for the recent star formation im- 
plied by the young ages. In an analysis of stellar absorption line gra- 
dients in the bulges of 32 nearby edge-on spirals (SO - Sc), Jablonka 
et al. (2007) also find that, when compared at a at a given veloc- 
ity dispersion, galaxy bulges resemble elliptical galaxies. Recently, 
MacArthur et al. (2008) have extended these results to bulges at 
intermediate redshifts (0. 1 < z < 1) finding that the mass assembly 
history of bulges, inferred via the Fundamental Plane, is indistin- 
guishable from that of pure spheroidal galaxies (E/SOs) at a given 
mass. 

A significant outlier from this trend, however, is the bulge of 
our own Milky Way (MW). Detailed photometric (Zoccali et al. 
2003) and spectroscopic (Rich & Origlia 2005; Zoccali et al. 2006) 
studies of abundance ratios of individual MW bulge stars reveal old 
and a-enhanced SPs that must have formed long ago and on short 
timescales, and no evidence is found for the presence of a younger 
SP. This observation is difficult to reconcile with the MW's small 
ao and the presence of a stellar bar (e.g. Lopez-Corredoira et al. 
2007). As Thomas & Davies point out, the MW studies sample a 
larger physical radius than those of external galaxies so this conun- 
drum could be resolved if there is a positive age gradient in the 
bulge. 

Adding further complexity to the picture, in a study of 
absorption-line maps of 24 early-type spirals, Peletier et al. (2007) 
find that while inclined samples, which sample the disk-free outer 
bulge regions, reveal uniformly old SPs, randomly oriented sam- 
ples, spanning the same range in Hubble type, show a wide range 
in SP parameters. This discrepancy, along with their observation 
of central dips in the velocity dispersion profiles of half of their 
bulges, indicative of the presence of a colder central component, 
can be reconciled in a picture whereby the centers of most early- 
type spirals contain multiple kinematic components: an old and 
slowly rotating elliptical-like component, and one or more disk- 
like, rotationally supported, components which are typically young, 
but can be old. The sum of all central components make up the pho- 
tometrically defined bulge, and the different components contribute 
differently to the integrated light depending on the viewing angle. 
Their conclusion of a rejuvenation of the central region through 
dynamical disk-like processes, however, is difficult to understand 
in the context of the above results in that the "rejuvenated" low-cr 
ellipticals do not have disk components. 

At the heart of these conflicts is the fact that extant studies of 
local bulges are limited in one way or another (sample size, range 
in Hubble type, data of insufficient depth and/or spatial resolution 
for studies of gradients and assessment of disk contamination, etc.). 
Our aim is to remedy this situation with a systematic and homoge- 
neous study of the SPs and kinematics of bulges for the full range 
of Hubble types. A first, and important, step is to establish a reliable 
procedure for extracting a proper representation of the SP content 
that does not suffer from the degeneracies and limitations discussed 
above. In this paper, we present a pilot sample of long-slit optical 
spectroscopy of eight nearby, low-inclination, spiral galaxies ob- 
served with Gemini/GMOS with which we develop such a tech- 
nique and attempt to address the formation of spiral bulges through 
a study of their radial profiles in SPs and kinematics. 

The organization of the paper is as follows. In g2] we detail the 
sample selection, observational strategy, and reduction procedures. 



Sj3] outlines our technique for extracting SP parameters from inte- 
grated spectra including details of the fitting algorithm, providing 
specifics about the weighting of individual pixels and our so-called 
"cr-clipping" procedure for masking deviant pixels, as well as a de- 
scription of the model SP templates used in the fits. Our results are 
presented in Sj4]which includes the derived SP parameters and their 
radial profiles, kinematic profiles (rotation velocity and velocity 
dispersion), an investigation of correlations of central parameters, 
and a comparison with previous studies. We discuss in Sj3] results 
for each individual galaxy in the context of bulge formation scenar- 
ios, and ^summarizes the salient points from the entire analysis. 
Finally, we provide several appendical sections: ijA] presents the 
observed radially resolved spectra for our galaxies. ^B] discusses 
derivations of SP parameters for our spiral galaxies based on the 
Lick-index system and of full spectrum fitting of single SSPs. All 
techniques are compared and contrasted in !jC] 



2 DATA 

The long-slit spectroscopic data for this study were collected using 
the Gemini Multi-Object Spectrograph (GMOS; Hook et al. 2004) 
on the 8-m Gemini North telescope at Mauna Kea in Hawaii. The 
GMOS detector consists of three 2048x4608 CCDs with 13.5 ^m 
pixels providing a spatial resolution of 0.072"/pix and a dispersion 
of 0.45 A/pix with the B600_G5303 grating. The internal stability 
of the slit mask requires two equally spaced bridges between the slit 
edges which lie along the spatial direction. The spacings between 
the 3 CCD detectors manifest as small (~17 A) gaps in wavelength. 

The high sensitivity of the B600_G5303 grating in the blue 
matches well our simultaneous spectral coverage of ~4050- 
6750 A. This range includes most of the major atomic and molecu- 
lar features to disentangle age and metallicity effects in integrated 
galaxy spectra (Worthey 1994). The slit field of view (FOV) was 
5' (length) x 2" (width) for all of our observations. The choice of 
slit-width was prescribed by the need to maximize signal-to-noise 
(S/N) in the outer disk while maintaining adequate spectral resolu- 
tion throughout. The 2" slit and B600 grating give a top-hat 10.8 1 A 
full width at half maximum (FWHM) resolution and a Gaussian in- 
strumental resolution of 0.8 ± 0.02 A (as measured from the width 
of the narrowest sky emission lines). 



2.1 Galaxy Sample 

Ideally, we desire a large and homogeneous sample of nearby 
spiral galaxies spanning the full Hubble sequence to study spectro- 
scopically the SPs of spiral bulges and disks. The sample should 
also include barred and non-barred spirals in order to test for the 
expected mixing effects by a bar. The size of the current sample 
was a compromise between a reasonable telescope time request for 
a pilot study and the need to measure systematic variations within 
that sample. We have thus narrowed in on the following criteria: 

• Hubble type Sa-Sd (with emphasis on later types) 

• Mix of barred/unbarred systems and SB profile types 

• Face-on (inclination < 35°) 

• Blue Galactic extinction A_g =4 x E(B — V) < 0.5 mag 
(Schlegel et al. 1998) 

• Blue major axis <2x slit length (< 10' 5), but large enough to 
resolve bulge and inner disk (> 2"). 

The galaxy sample and catalog information from the NASA 
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Table 1. Galaxy Sample: catalog information. 



Names Hubble Diam RA DEC Vhelio M s A s 

NGC/1C UGC Type a(') X b(') (J2000) (tar 1 ) (mag) (mag) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



N0173 


U369 


SA(rs)c 


3.2 


X 


2.6 


00 h 37 m 12; 


'47 


+01°56'32'.'l 


4366 


13.70 


0.110 


N0628 


U1149 


SA(s)c 


10.5 


X 


9.5 


01 h 36 m 41 ! 


'77 


+15°47'00'.'5 


657 


9.95 


0.301 


N1015 


U2124 


SB(r)a 


2.6 


X 


2.6 


02 h 38 m ll ? 56 


-01 o 19'07'.'3 


2629 


12.98 


0.139 


N7490 


U12379 


Sbc 


2.8 


X 


2.6 


23 h 07 ra 25 1 


' 17 


+32°22'30 / /2 


6213 


13.05 


0.362 


N7495 


U12391 


SAB(s)c 


1.8 


X 


1.7 


23 h 08 m 57! 


18 


+12°02'52"9 


4887 


13.73 


0.371 


N7610 


U12511 


Scd 


2.5 


X 


1.9 


23 h 19 m 41 ! 


'37 


+10°11'06'.'0 


3554 


13.44 


0.171 


N7741 


U12754 


SB(s)cd 


4.4 


X 


3.0 


23 h 43 m 54> 


'37 


+26°04'32{'2 


751 


11.84 


0.323 


10239 


U2080 


SAB(rs)cd 


4.6 


X 


4.2 


02 h 36 m 27 t 


'88 


+38°58'liy7 


903 


11.80 


0.307 



Notes — Cols. (1) & (2): NGC/IC, and UGC Galaxy IDs. Names highlighted in bold are those used subsequently to refer to the particular galaxy. Col. (3): 
Hubble Type (from NED). Cols. (4 & 5): Major (a) and minor (b) axis diameters (in arcmin). Cols. (6 & 7): Galaxy RA and DEC (J2000). Col. (8): Heliocentric 
radial velocity. Col. (9): Total 73-band magnitude (Vega) from NED. Col. (10): B-band galactic extinction from Schlegel, Finkbeiner, & Davis (1998). 



Extragalactic Database (ned£| are shown in Table[T| Optical imag- 
ing for all eight galaxies was collected using the Large Format 
Camera on the 200-inch Hale telescope at the Palomar Observatory 
in December 2006. Observations were obtained in the BVRI Bessel 
filters for all galaxies except N0173, for which only B&V were ob- 
tained (due to bad weather). Surface brightness profiles were ex- 
tracted from these images and bulge and disk photometric parame- 
ters were derived from bulge-to-disk (B/D) decompositions into an 
exponential disk plus a generalized power-law (Sersic) profile for 
the bulge, as listed in Table [2] We refer the reader to MacArthur, 
Courteau, & Holtzman (2003, hereafter Mac03) for details about 
profile extraction and B/D decompositions. Note that here we refer 
to the disk scale length as rd and the bulge effective (half-light) 
radius as r e for a profile characterized by the Sersic n shape pa- 
rameter. 

As is typical for late-type bulges (e.g. de Jong 1996; Graham 
2001; Mac03), most of our decompositions reveal bulge profile 
shapes close to exponential, which translates to a Sersic index n = 1. 
Most notably, the bulge of N7495 was best fit with n = 2.3. Bulges 
with n > 2 are typically found in earlier-type bulges thought to be 
more akin to elliptical galaxies (e.g. Fisher & Drory 2008). This 
suggests perhaps that the Hubble type of this galaxy is earlier than 
Sc. In fact, as noted in NED, the original UGC catalogue does refer 
to this as an "early spiral" (Nilson 1973). On the other hand, our 
only Sa bulge, U2124, is best fit with n = 1.5 which is low for an 
early-type bulge. However, the strong bar in this galaxy renders the 
decomposition less certain. 



2.2 Standard Star Sample 

We observed seven Lick/IDS standard stars (Worthey et al. 1994) 
with the same observational set-up as the main targets for flux 
calibration of the galaxy spectra. These included dwarfs and 
giants with a range of physical parameters relevant for star- 
forming galaxies (spectral type: A1-F5; effective temperature: 
T ef f ~ 4700-9000 K, surface gravity: log(^) = 2.4^1.3; and metal- 
licity: [Fe/H] = —0.45-0.07). Flux calibration was achieved by 
comparison of our standard star observations with flux calibrated 
spectra of stars in common with the ELODIE archive (Moultaka 
et al. 2004; see http://atlas.obs-hp.fr/elodie). 



1 http://nedwww.ipac.caltech.edu/ 



Table 2. Bulge and Disk Photometric 
Scale Parameters 



Name 


r e 


Sersic 




NGC 


(") 


n 


(") 


(1) 


(2) 


(3) 


(4) 


N0173 


3.8 


1.5 


21.9 


N0628 


10.5 


1.1 


68.2 


U2124 


5.7 


1.5 


12.2 


N7490 


4.5 


1.6 


21.4 


N7495 


3.6 


2.3 


17.2 


N7610 


1.7 


0.8 


17.6 


N7741 


6.7 


1.0 


42.3 


10239 


8.2 


0.9 


43.8 



Notes — Photometric bulge and disk scale parameters 
based on optical imaging (either R or /-band, but 
V for N0173) from the Palomar 200-inch telescope. 
Errors on the disk scale lengths are < 10% and < 20% 
for bulge parameters. See Mac03 for details of our 
B/D decomposition scheme. Col. (1): Galaxy ID. Cols. 
(2-4): B/D decomposition values for bulge effective 
radius (r e ), Sersic n, and disk scale length (r d ). 

2.3 Integration Times and Observing Strategy 

The on-target integrations were divided into separate exposures to 
filter out cosmic rays and confirm weak features. On-target In- 
tegration times were 3 x 15 min for Mb < 10.5, 3x18 min for 
10.5 < A1 B < 12, and 4x18 min for M B > 12, for total integra- 
tions of 45-72 min per galaxy. An estimate of the sky background 
was measured from the slit edges for galaxies with major axis di- 
ameters < 4'. For galaxies wider than the slit, we divided our expo- 
sures into sky/galaxy/sky sequences with integration times for sky 
(off-target) pointings of ~l/5 of the on-target integration^. See 
Fig. [T] for the observational set-up for each galaxy. Sky positions, 
when necessary, are shown by the yellow slits. Sky exposures were 
binned to achieve the same S/N as the galaxy. 

2.4 Data Reduction 

The relevant measurements for this study do not require absolute 
spectrophotometric calibration of the data. However, they do re- 

2 In hindsight, sky integration times of 1/3-1/2 of the on-target times would 
have been preferred. The short sky exposures prevented a deeper explo- 
ration of the low-SB disk regions. 
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Figure 1. Observational set-up for the eight galaxies in our sample. The background images are from the Canadian Astronomy Data Centre's Digitized Sky 
Server (CADC; http://cadcwww.dao.nrc.ca/ 1. The blue line represents the slit, the red (dashed) box and long arm represent the FOV of the GMOS wavefront 
sensor camera, with the box at the end of the arm centered on the guide star. The panels for large galaxies (NGC 628, NGC 7741, and IC 239) also show, as 
yellow lines, the sky offset positions. The FOV for the CADC pictures differ in all the panels but the slit length is everywhere the same (5'). 



quire extreme care in removing instrumental effects in order to 
achieve accurate relative calibration along both the dispersion and 
spatial directions. Proper subtraction of the sky and instrumental 
and far-field scattered light is also crucial. Finally, accurate flux 
calibration is required for the galaxy observations, particularly for 
long-baseline measurements such as those of molecular indices and 
full spectrum fitting. We went to great lengths to assure the highest 
possible quality in our final data products. In the interest of space, 
the detailed data reduction procedures will be described in a forth- 
coming paper. 



3 EXTRACTING STELLAR POPULATION 
PARAMETERS FROM SPECTRA 

Many different techniques can be used to characterize the stellar 
population parameters of an integrated (optical) spectrum. We have 
explored several of these in detail. Appendix ijBjpresents our anal- 
yses based on the Lick system of absorption-line indices and that 
of full spectrum single SSP fits. Both methods produce SP param- 
eters based on SSPs and are thus extremely sensitive to the most 
recent episode of SF, even when its contribution to the SP in terms 
of mass is insignificant. The line-index technique is also very sen- 
sitive to any ongoing SF as the emission-lines from H II regions fill 
in the underlying absorption features, thus thwarting reliable age 
estimates. We conclude in Appendix SjC]that the SPs of our spiral 
galaxies are not well represented by SSPs and an accurate represen- 
tation must represent true average values of the SP parameters, i.e. 
integrated over the SFH of the galaxy. In the following sections we 
describe our devised method for achieving this goal, which refer to 
as "full population synthesis". 



3.1 Full Population Synthesis 

Population synthesis is a technique applied to the integrated light 
of composite systems to determine the fraction of light (or mass) 
contributed by each of the main stellar groups comprising the sys- 
tem. Using a linear combination of stellar or model templates that 
cover the optical spectral range, the temperature mix of a galax- 
ian (composite) system can be determined and used to explore its 
evolutionary history and to predict the flux of the system at other 
wavelengths as a test of the synthesis (Jacoby, Hunter, & Chris- 
tian 1984; Pickels 1985; Heavens et al. 2000; Cid Fernandes 2007; 
Tojeiro et al. 2007). 

The "full population synthesis" technique developed here con- 
sists of deconvolving an observed galaxy spectrum into basic com- 
ponents representing the relative contributions of a basis set of tem- 
plate spectra. It is important for the method of decomposition to 
allow for constraints to be placed on the computed parameters for 
which we turn to the method of bound constrained optimization 
(Byrd et al. 1995; Zhu, Byrd, & Nocedal 1994, 1997). The problem 
can be summarized as solving for 

min/(a;) 
subject to I ^ x ^ u ' 

where f(x) is the merit function which depends non-linearly, in 
general, on the fit parameters Xj, and lj and Uj are lower and upper 
bounds respectively on each parameter Xj . Additionally, the gradi- 
ent g of the merit function / with respect to the fit parameters must 
be available. The algorithm performs the minimization iteratively 
where, at each iteration, the current fit parameters are varied by 
amounts determined by the current partial derivatives of the merit 
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function. To start the iteration an initial configuration must be sup- 
plied. 

Ideally, the locus of the merit function's minimum will be a 
single point in the vector space defined by the basis library of tem- 
plate spectra. The necessary conditions for this ideal case are: (1) 
the basis set is complete and non-redundant and (2) there is no noise 
in the data. In this situation the solution will be independent of the 
initial configuration. However, in any real synthesis we are limited 
by observational uncertainties, resolution, physical differences be- 
tween the observed spectrum and model templates (e.g. abundance 
ratios, dust, emission lines from H II regions, any non-stellar con- 
tribution from AGN, etc.) and possible degeneracies in the library 
itself. Hence, the locus of the minimum merit function will be a fi- 
nite volume in the basis space. The ultimate success of a population 
synthesis is therefore limited by the observational accuracy of the 
spectra and the match between data and template library. 



3.2 The Program 

The problem of bound constrained optimization arises in numerous 
scientific and technical applications and has been an active area of 
research in the field of numerical analysis. As a result, there are 
now a number of robust optimization routines available; XRQP 
(Bartholomew-Biggs 1979) and LANCELOT (Conn, Gould, & 
Toint 1992), to name a few. The program of choice for the present 
problem is that of Zhu et al. (1994; 1997) called L-BFGS-E0. This 
is a FORTRAN 77 code that solves large non-linear optimization 
problems with simple bounds. It has been well tested and is ideal 
for the current problem. A full description of the L-BFGS-B algo- 
rithm is beyond the scope of this paper, but the interested reader is 
referred to Byrd et al. (1995) and Zhu et al. (1994; 1997). 

A driver program popsynth was written which calls the L- 
BFGS-B routine, popsynth performs the following steps: 

(i) Set up input values for call to setulb, the calling statement 
of L-BFGS-B 

(ii) Set up initial values for fit parameters Xj 

(iii) Read in the galaxy spectra and associated variance vectors, 
and the template library 

(iv) Compute the merit function and its gradient for the given 
set of fit parameters 

(v) Call the setulb subroutine with the new merit function and 
gradient values 

(vi) Test for convergence, if not repeat (iv)-(vi) 

(vii) Output stellar fractions to data files. 

Zhu et al. (1997) provide a complete description of the input param- 
eters, but a brief introduction is relevant here in order to understand 
the convergence criteria of the algorithm. 

The input variable factr is a double precision number 
which defines a tolerance level in a termination test for the algo- 
rithm. The iteration will stop when 



fk — fk-l 



max(\f k \, |/j 



k + l | 



^ factr * epsmch, 



where epsmch is the machine precision which is automatically 
generated by the L-BFGS-B code (equal to 2.22 x 10~ 16 on the 
Linux machine used for this project.) A value of f actr = 10 7 was 



chosen for "moderate" accuracy. Similarly, pgtol is a tolerance 
level for termination of an iteration if 

maa;(jproj gi\, i = l..n) ^ pgtol 

where proj gt is the i th component of the projected gradient. For 
the current application pgtol was set to 10~ 5 . The iteration ter- 
minated according to the latter criterion for all tests and runs of the 
program. 

The merit function to be minimized is defined as 



Merit 



N — M 



where 

Ri = Wi(Gi — Si 



(2) 



(3) 



In the above expressions, N is the number of data points included 
in the fit, M is the number of templates in the library, Gi is the 
observed galaxy flux at wavelength Xi, Wi is the weight of the i th 
pixel (see ^3.3b . and Si is the modeled galaxy flux at Xi given by 



(4) 



3=1 



3 Can be downloaded from http://www.ece.northwestern.edu/'^nocedal/lbfgsb. 
as of writing. 



where Fji denotes the flux at the i th wavelength of the j th tem- 
plate, and Xj is the relative contribution of the j template to the 
synthesized spectrum. 

In order to prevent negative stellar contributions, as a non- 
physical artifact of the decomposition, the spectral coefficients are 
subject to a lower bound on the fit parameters such that 



3.3 Weighting 

The nature of our fitting algorithm is particularly useful for mask- 
ing out undesirable regions in the spectrum. The following patholo- 
gies, which differ for each galaxy due, in most part, to velocity dif- 
ferences, are examples of regions that one may want to omit: (i) 
gaps in wavelength coverage due to the spacing between the three 
GMOS CCDs, (ii) strong and variable night sky-line regions, and 
(iii) emission line regions (due to H II regions and/or planetary neb- 
ulae). As a first guess, all of (i), (ii), & (iii) are masked for the initial 
fit of each spectrum. Subsequently, we allow for all points to enter 
back into the fit (subject to the constraints defined below) except for 
the gap regions, in which no real data exist, thus are always fixed 
to zero weight. 

3.3.1 a Clipping 

We define a mask value, rrn, for each datum such that 
( 1 , include in current fit 

771 ' — s 

1 \ 0, remove from current fit 
and the weight, Wi, as, 



(5) 



where Oi is the error on Gi. 

In a general weighted fit, very dissimilar weights can be dis- 
tributed among the data points. As a result, the number of datum, 
htmlV, does not always relate directly to the actual information feed- 



ing the fit. Instead, an effective number of points, N e f / , can be 
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estimated by the sum of the weights divided by the mean weight 
(Gonzalez 1993) (note that in the following, all sums run from 
i = 1..JV, where N is the total number of data points): 



N, 



eff = 



(£'• 



K> (E 1 



0/E 1 



E'< 



(6) 



Next we define the effective number of templates entering 
each fit, M e ff . This is required as it is possible for a significant 
fraction of the total number of templates considered to end up with 
zero contribution to the fit. M e f / is defined as follows (here sums 
run from j = 1..M, where M is the total number of templates con- 
sidered): 



eff 



E*i 



E' : ' E ' . E ' E' 



where Fj is the fractional contribution of SSP model Mj 



(7) 
(thus 

^2 Fj = 1 by definition). 

The expected variance (of the mean) based on the measure 
ment errors is 

2 N B ff 



E w « 

and the typical variance from the fit is, 

E( w * A ?) 



2 



AT, 



AT, 



M, 



e// 



E< 



(8) 



(9) 



where A* = |G< - 5»|. 

Finally, the fit ^ 2 is the ratio of the mean typical and expected 
variances, 



2 _ 2 / 2 
X — °~fitl°~exf> 



(10) 



At each iteration, the quantities in Eqs.[4}j9]are computed and 
the individual masks, mj, are assigned according to the criterion, 

= f 1, if Aj < AT CT o- /rt 
1 1 0, otherwise 

On the first iteration we set 7V CT = 5, the second has N a = 3.2, and 
two final iterations have N a = 2.8. We refer to this iterative masking 
scheme as "cr-clipping". 

For our data, the typical number of data points for each spec- 
trum is in the range N = 1280-1340; after the "cr-clipping" itera- 
tions described above, the typical effective number of points in the 
final fit spans the range N e /f = 850-1250. In this analysis we use 
M = 70 model templates (see £j3.4b and typical values for M e f / are 
in the range 2-6. 



3.4 Stellar Population Templates 

The stellar population templates used throughout this analysis are 
from the models of BC03 which provide spectra of simple stel- 
lar populations (SSPs), single-Z/single-bursts of star formation. 
These cover the wavelength range 3200-9500 A, a large range of 
ages (0.0001-20 Gyr), and metallicities (Z = 0.0001-0.05) at a re- 
ported resolution of ~ 3 AFWHM. We use here the models with 
the Chabrier (2003) initial mass function (IMF). 



3.4.1 Template Characterization 

Given the importance of knowing the precise resolution of the tem- 
plates for comparison with galaxy data and for measuring the small 
velocity dispersions of spiral galaxies (see we undertook a 

deep investigation into the true resolution of the BC03 models, 



which use the spectral library STELIB of Le Borgne et al. (2003), 
for a range of model SSPs and allowing for a wavelength dependent 
(velocity-like, i.e. constant with A A/A) term in addition to a term 
that is constant with wavelength. Our analysis is carried out using 
an upgraded version of the Mo ve 1 algorithm described in Gonzalez 
(1993). In brief, Movel derives velocity shift, velocity dispersion, 
and 7 (a measure of the relative line strengths, with 7=1 indicat- 
ing a perfect match) through a Fourier transform procedure. In ve- 
locity space, the observed galaxy spectrum can be characterized 
by the spectrum of its population convolved with the instrumental 
response function and a broadening function which characterizes 
the distribution of stellar radial velocities along the line of sight 
(LOS VD). The "population" is represented by a template (observed 
or model) of higher resolution than the galaxy spectrum. Movel 
derives the broadening function that best matches the template to 
the galaxy spectrum. To then infer the absolute rotation and veloc- 
ity dispersion parameters of the observed spectrum, one must know 
very accurately the characterization of the template used. 

For the purpose of characterizing the BC03 resolution (both 
globally and as a function of wavelength), we use as a template 
spectrum the very high-resolution (effectively infinite for this ap- 
plication) solar spectrum of Wallace, Hinkle, & Livingston (1998) 
along with the BC03 model SSPs, which are each effectively dif- 
ferent coadditions (in relative weight and stars used) of stars in 
STELIB. The solar template and model SSPs are processed through 
Movel in 200 A wavelength intervals centered around 26 spectral 
features, from 3575 A up to 7400 A defined by 21 Lick-like indices 
plus the Ha and Ca II H&K lines, and 3 more somewhat ad hoc re- 
gions shortward of Ca II H&K and longward of Ha to cover the full 
optical range. Spectral features in any given SSP that are not well 
characterized by the solar template are weighted out accordingly. 
We then fit the trend of resolution vs. A for a large range of SSP 
templates using three different constraints; (i) fitting for constant 
terms with both A and velocity with both terms as free parameters, 
(ii) fixing FWHM = 3 A and fitting for a velocity term, and (iii) forc- 
ing only a linear FWHM component. This procedure revealed that 
the resolution trend of the BC03 models is much better represented 
by a FWHM constant with A plus a velocity term (constant with 
A A/ A). 

To investigate the cause of this discrepancy, we downloaded 
the STELIB librarjQ and examined the individual spectra. First we 
noted that the starting A and dispersion are not all 3200 & 1 A as 
stated in the documentation on their website. The reason is that 
the 2003 v3.2 library stars had their radial velocities removed, but 
the spectra were not rebinned to the Ao = 3200 & disp = 1 A. A 
look into the headers revealed a spread in the starting wavelength 
of RMS Ao = 1.09 A and in dispersion of RMS disp = 0.00033 A. As 
these two differences in the relative wavelength calibration are cor- 
related and, if fixed values are assumed when coadding the indi- 
vidual stars, the combined effect will degrade the effective reso- 
lution of the model SSPs. In measuring the resolution as a func- 
tion of wavelength for 12 individual solar-type stars in the STELIB 
library we found, not surprisingly, that the individual stars are 
consistent with a pure constant resolution (no velocity-like broad- 
ening at all) and we measured the best constant FWHM to be 
3.08 A(± 0.03). We thus concluded that the above-mentioned Ao 
& dispersion mismatch was not taken into account when the stars 
were coadded to create the BC03 model SSPs, thus resulting in a 
velocity-like broadening in addition to the constant resolution term. 



http://www.ast.obs-mip.fr/users/leborgne/stelib/index.html 
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We also found a spread in the velocities of the 12 stars (whose 
measured relative velocities were removed in STELIB) with an 
RMS= 11.4kms -1 , which will also contribute to the velocity-like 
broadening of the SSPs. 

The above procedure leads to the following conclusions about 
the true effective resolution of the BC03 SSPs: (i) The actual spec- 
tral resolution of the models is worse than 3.0 A (significantly 
larger than for individual stars in STELIB); (ii) The best con- 
stant (no velocity broadening) FWHM is 3.40 A(± 0.04); (iii) If 
the 3.08 A resolution derived for individual stars is assumed, the 
SSP models (co-added STELIB stars) need a significant additional 
velocity broadening (35.8 ±3.6 km s~ ), but the fit is not as good 
as when assuming a constant ~3.4A resolution. After exploring 
the fit results under all the above assumptions, the adopted (best) 
values for the FWHM resolution and velocity broadening for the 
BC03 models are found to be: 

FWHMbcos = 3.375 A & <r B C03 = 11.40 kms -1 . 

These results can be understood as the combination (convolution) 
of: (i) the 3.08 A FWHM resolution of each individual star, plus (ii) 
the 1 1 .40 km s ~ 1 radial velocity RMS error in STELIB , and (iii) an 
additional 1.38 A residual FWHM degradation in the models due 
to the errors in the assumed initial wavelength and dispersions in 
STELIB in the coaddition of BC03. 

Finally, we used the same solar template to map the 
wavelength-scale of STELIB. The systematic deviations are no 
larger than AA = — 0.5 to +0.5 A across the observed wavelength 
range, yet are relevant enough for a fair comparison with high- 
quality modern data. These true systematic trends in resolution 
and wavelength scale of the BC03/STELIB spectra are in excel- 
lent agreement with those found in the similar analysis by Koleva 
et al. (2008). 

In our full population synthesis fitting, all of the above is ac- 
counted for. In particular, we first distort our observed galaxy wave- 
length scale to match that of STELIB. We then account for resolu- 
tion effects by convolving all model and galaxy spectra (at all radii) 
to the same resolution: top-hat = 10.81 A, a=168kms _1 . These 
account for the slit-width and approximately highest velocity dis- 
persion in our galaxy sample (measured on the first iteration and 
limited to the central, high S/N regions), respectively. This disper- 
sion is sufficiently greater than that of the native BC03 models such 
that any residual variations, e.g. due to different combinations of 
stars in different SSPs, will be smoothed out. We discuss in i ]4.4l our 
accounting for the BC03 effective resolution for kinematic mea- 
surements, which are iteratively derived from the full spatial and 
spectral resolution data. 



3.4.2 SSP Template Library 

We are limited here to a template library that comprises a discrete 
sampling of SSP ages and metallicities. Ultimately, we would like 
our library to be a complete and non-redundant basis set, but this 
would require an exact match between models and data that is not 
realized in practice. Through a series of tests and careful examina- 
tion of the individual SSP spectra, we converged on a selection of 
SSPs at 14 ages, roughly logarithmically spacecjf], and 5 metallic- 
ities (excluding only the lowest metallicity provided in the BC03 



5 The spacing in not exactly logarithmic, particularly at ages of about 
0.004 Gyr at the lowest Zs, where we noted anomalous SSP spectra in the 
BC03 models. The spectra around this age were inspected individually to 



models which is the least well modeled and exhibits anomalous 
SSP spectra). This provides a total basis set of 70 SSP spectra, all 
shown in Fie.ljFi 

We are interested in both the light- and mass-weighted con- 
tributions of each SSP to the integrated galaxian light. All SSPs 
in our library have been normalized to their y-band flux, so the 
V-band light-weight for SSPj is computed as (here sums run from 
j = 1 ..M, where M = 70 is the total number of SSP templates con- 
sidered): 



il 



i 



and the mass-weight as 



= x j *(M/L v ) j /y(x j *(M/L v ) j ) 



3 



(ID 



(12) 



where (M /L v )j is the V-band mass-to-light ratio of the j th SSP. 
Thus the average light- weighted quantities are computed as: 



(A), A* 



and 



(Z)i = ^2h* Zj 



(13) 



and the mass-weighted as: 

(A) m = m,j * Aj and (Z) m = 



2^m j *Z j . (14) 



Finally, we also compute the stellar I/-band M/L ratio of the com- 
posite system as M/Ly =^2, (h * (M/Lv)j). 

A number of caveats exist when applying model SSPs to the 
integrated light of a star forming galaxy. First, the models are based 
on empirical libraries of stars which are necessarily limited to the 
range of T e //, [Fe/H], log(<?), and [a/Fe] of stars within the solar 
neighborhood. Stellar populations in other regions in our Galaxy 
and in external galaxies will not be confined to this parameter space 
and thus not be faithfully represented. However, extrapolations in 
the (T e ff, [Fe/H], log(g))-space within the models do provide ex- 
tensive coverage, with the extremes (low T e //, small [Fe/H]) be- 
ing the least well modeled. Additionally, the BC03 models do not 
include the effects of non-solar abundance ratios and a full assess- 
ment of potential biases awaits the upcoming implementations of 
stellar population models which include abundance pattern varia- 
tions. Second, the models do not account for emission from plan- 
etary nebulae and H II regions which can be dominant in actively 
star forming galaxies such as the late-type spirals in our sample. In- 
deed, many of our galaxy spectra show strong emission lines (e.g. 
H/3, [O III] AA5007,4959, Ha, [N II], [S II], etc.; see Fig.|AT). We 
successfully circumvent this issue with our iterative masking pro- 
cedure described in ^13.3.11 Third, and again particularly important 
in star forming systems, the optical stellar light will likely suffer 
from the extinction and reddening effects of dust. While there is 
no straightforward way to account for dust effects on integrated 
spectra, we attempt to determine if a particular spectrum is affected 
by extinction by comparing dust-free fits with those incorporat- 
ing a model for dust extinction. We adopt the two-component dust 
model of Chariot & Fall (2000, hereafter CF00). The CF00 model 
accounts for the finite lifetime (~10 7 yr) of stellar birth clouds 



find the "best behaved" age (i.e. bluer for younger ages). The 0.006 Gyr 
SSPs had the most reasonable spectra. 

6 Note that we are not concerned here with model age predictions that are 
older than the age of the Universe. Model ages are not precisely calibrated 
and, as such, we are concerned primarily with relative trends. 
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Figure 2. Spectra of the 70 SSP templates used in the population synthesis fits. Metallicity increases from left to right. Different ages have different colours 
and are labeled in leftmost panel. All spectra are normalized to their V-band flux. 



through the parameter fi, which represents the fraction of total dust 
absorption contributed by diffuse interstellar medium (cirrus) dust. 
The wavelength dependence of the effective absorption curve is 
oc A -0 ' 7 , so the "effective" extinction in a given band, f\, is given 

by, 



f tv(A/5500 A)" ' 7 for t s? 10 7 yr, 
\ //fy (A/5500 A)" ' 7 for t > 10 7 yr, 



where t is the age of any single stellar generation. 

For each SSP spectrum in our library, we create extincted ver- 
sions with a fixed /i = 0.3 (the average value for star forming galax- 
ies) and varying fv in steps of 0.5. For reference, in this model, 
tv = 2 is equivalent to a colour excess E(B — V) ~ 0.1 for a so- 
lar metallicity template (see Fig. 1 in MacArthur 2005). We per- 
form the full synthesis first on the dust-free SSPs. The results are 
recorded and the Merit function (Eq. [2} tabulated. The fits are 
then performed on the extincted SSPs, starting with fy = 0.5. If 
the Merit function has decreased from the previous fit, the next 
extinction level is fit, until a minimum in the Merit function is 
reached. While the reddening effects of dust add an extra degen- 
eracy to the fits, we always reach a convergence (minimum Merit 
in the Age/Z/dust-space). 



4 RESULTS 

4.1 Full Population Synthesis Fits 

We here present example fits from the code and procedure de- 
scribed in § 13,11 Fig.[3]shows the central, r = 0", spectrum fit for all 
eight galaxies. In each plot, black and red lines are respectively the 
data and synthesis fits. The gray shading indicates regions masked 
in the fits as determined by our "cr-clipping" procedure described 
in § 13.3.11 Shown at lower right on each panel are the average light- 
weighted age, (A); , and metallicity, (Z); , effective fv, and x 2 of 
the fit. The bottom panels show the percent data— model residuals. 
In all cases, the full synthesis fits represent the galaxy spectra re- 
markably well. Our cr-clipping procedure clearly does a good job 
at masking out regions where emission lines are present. 

Closer examination of the data— model residuals for all of the 
spectra that show little or no evidence for current SF reveal the 
same overall rolling shape, with the most significant feature occur- 
ring around A ~ 5150-5450, a region rich with Fe and Mg absorp- 
tion lines. While this may be due to a flux calibration problem, it 
could also be an indication of non-solar abundance ratios that are 
not taken into account in the BC03 models. A similar pattern was 
observed in the fit residuals of SDSS galaxy spectra in the studies 
of Panter et al. (2007) and Cid Fernandes (2007). The fact that these 
independent data sets show the same features effectively rules out 
a flux calibration problem. Panter et al. (2007) suggest that these 
features may be due to an error in the balance of K-M giants in 
the BC03 models, whereas Cid Femandes (2007) attributes them to 
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Figure 3. Central observed spectra (black) and full population synthesis fit (red) for all eight galaxies. Gray shading indicates regions masked in the fit as 
determined by our iterative "cr-clipping" procedure (see £|3.3. It as well as the CCD gap regions (green vertical dash-dotted lines) which are always masked. 
Shown at lower right on each panel are the average light- weighted age, (A); , and metallicity, (Z); , effective ry, and x 2 of the fit. The bottom panels 
show the percent data— model residuals. In the bottom panels, dashed vertical lines indicate variable sky-lines (which can be difficult to account for during 
sky-subtraction), and dotted vertical lines indicate emission lines prevalent in star forming (H II) regions. 
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non-solar a-element abundance ratios. We may be able to test this 
in future work with new SSP models that include non-solar a/Fe. 

For now, however, the best diagnostic we have of a/Fe abun- 
dance comes from the Lick index measurements. In particular, in 
a plot of an index primarily sensitive to Fe versus one sensitive 
to an Qf-element, an enhancement in the latter can be identified 
if the data extend beyond the solar a/Fe model grids. We exam- 
ine possible a-enhancement in our data in Fig. [4] where we plot 
(Fe) vs. Mgb. For clarity, only the central, r — 0", data point is 
plotted for each galaxy (for direct comparison with central spec- 
tra shown in Fig. [3}- The black squares are measured directly from 
the galaxy spectra, whereas the pink squares are measured from the 
full synthesis model fits (red spectra in Fig. O. By definition, the 
pink squares must lie within the model grids, thus a difference be- 
tween the black and pink squares indicates that the models failed 
to reproduce the data simultaneously in this combination of nar- 
row index regions. Indeed a few galaxies show signs of super-solar 
a/Fe ratios in their centers, with the most significant enhancement 
in U2124, which also shows the most prominent features in its full 
synthesis fit residuals (Fig. [3}- This lends credence to the non-solar 
abundance ratio interpretation for the cause of the most prominent 
residuals that are not attributable to non-stellar radiation sources. 

Another possibility to examine the conjecture that the dom- 
inant features of the fit residuals in Fig. [5] are due to non-solar 
[a/Fe] uses the synthetic spectral library of Coelho et al. (2005) 
who present models for a range of T e //> log(g), and metallicity, 
here expressed as [Fe/H], with varying abundance ratios. The ef- 
fects of super-solar a/Fe are shown in Fig. [5] as the difference be- 
tween the solar and enhanced a/Fe models of Coelho et al. (2005) 
for T eff = 4000 K (top) and 5000 K (bottom), for a range of \og(g ) 
and [Fe/H]. The effects depend strongly on T e ff & log(g), but 
some of the features and wiggles are similar to those seen in our 
galaxy-model fit residuals, thus strengthening the conclusion that 
the wiggles are due, at least in part, to abundance ratio effects. Indi- 
vidual features such as those seen in Fig.|5]at A ~ 4150-4400, 4600, 
5 1 50-5450 (with a prominent feature at ~ 5 1 80 A), and 6200-6300 
A, are most evident in the residuals for the central spectrum of 
U2124 in Fig. [3] The issue of abundance ratio variations will be 
examined in more detail in a future work. 

Cid Fernandes (2007) also noted another striking feature in 
the data— model residuals of their mean spectrum for star-forming 
galaxies in the SDSS, selected on the basis of emission line 
diagnostic-diagrams. This feature takes the form of a shallow but 
broad "absorption" band around H/3, that appears whenever stars of 
~ 100 Myr contribute significantly to the spectrum. While Cid Fer- 
nandes points out that no straightforward physical explanation nor 
element identification for this trough exists, he also speculates that 
it is a side-effect of flux calibrations issues in the STELIB spectral 
library. We do not see such absorption troughs around H/9 for our 
star-forming galaxies (see, e.g., N7495 in Fig. [3}. We thus specu- 
late that the trough seen by Cid Fernandes could be a residual trend 
from the resolution issues of the BC03 models described in §3.4\ 
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Figure 4. (Fe) vs. Mgb as an indication of non-solar abundance ratios. Solid 
lines are iso-Z tracks with colour scheme: Z = 0.05 (red), 0.02 (Zq) (yel- 
low), 0.008 (green), 0.004 (light blue), 0.0004 (blue), & 0.0001 (purple). 
Gray-shaded dashed lines are iso-age tracks at ages of 1, 2, 5, 10, 13, & 
20Gyr (from light to dark shades; roughly left to right). The black squares 
are measured from the data. The pink squares are measured on the model fit 
spectra. 



the cumulative statistics of the contribution of each individual SSP 
template to the fits. Clearly, certain SSPs are favored, with many 
never entering any fit. This is partly the reason for the stability in 
the fits. With 70 templates to choose from, one might expect the 
fits to be very sensitive to initial conditions, for example. We have 
tested the fits using various initial conditions and the results are 
extremely consistent. Given that typically fewer than 10 templates 
have non-zero weight in the fit, the stability is not as surprising. 
Interestingly, the solar Z = 0.02 SSP rarely makes it into any of the 
fits and only at the oldest three age bins, whereas all other Zs have 
some contribution at a range of ages. While it is difficult to identify 
the exact reason for this (e.g. perhaps the real parameter space is 
not fully represented in the models, or some systematic in the data 
forces it to certain areas of the model parameter space), an investi- 
gation based on fits to other SPS models may provide some insight, 
and will be the subject of a future investigation. 



4.1.1 Relative Contribution of Individual Templates 

It is interesting to investigate which of the 70 SSP templates are 
most represented in the fits. We have no reason a priori to think that 
the SSPs at any given age and Z will be more likely to represent the 
SPs of our spiral galaxy data, but it is possible that certain model 
SSPs suffer some systematic anomaly rendering them poor repre- 
sentations of real SPs (we certainly observed some such patholo- 
gies when selecting our 70 template SSPs, see £13 ,4b . Table[3]shows 



4.1.2 Reliability of the Fits 

To test the robustness of our method and results, and to associate 
reliable error estimates to the derived individual age, Z, and ry 
estimates, we perform a Monte Carlo (MC) analysis on a subset 
of radial bins for each galaxy. A full-blown MC analysis for every 
spectrum in our sample is computationally prohibitive, so the sub- 
set analyzed was selected to be representative of the entire data set 
spanning the full range of S/N/A, SED shape, emission strength, 
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Table 3. Statistics of Individual SSP Templates Included in Full Population Synthesis Fits of Integrated Galaxy Spectra 
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Figure 5. Comparison of solar vs. [a/Fe] = +0.4 from models of Coelho 
et al. (2005), plotted as the %-difference of enhanced minus solar-scaled 
models, for effective temperatures of 4000 K (top) and 5000 K (bottom), and 
for a range of surface gravity and metallicity (as labeled in bottom panel). 
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Figure 6. Results from Monte Carlo simulations. Plotted are the logarithmic 
differences between the parameters measured in each MC realization and 
the best-fit value. Point types and colours distinguish the different galaxies 
(see legend of Fig. lCU . 



and radial coverage of each galaxy. Twenty realizations for each 
selected spectrum, after adding random noise, were then processed 
through our fitting algorithm in an identical fashion as the origi- 
nal galaxy data. Characteristic errors for the derived parameters are 
then taken as half the interval containing 68% of the MC realiza- 
tion results (i.e., the 1-er confidence interval). Fig. [6] summarizes 
the results from this exercise, showing the difference between the 



parameters derived from each MC realization relative to the best- 
fit value. There is evidence for a mild age/Z error correlation, and 
where not present, it turns up as an error degeneracy with fy. The 
errors are typically << 10% and, as expected, tend to increase with 
decreasing age, Z, and S/N/A. Error bars from these simulations 
are plotted in all subsequent figures. In particular, Fig.[8]reveals the 
radial sampling of the MC analysis. 

Beyond the small errors of the mean SP parameters, the MC 
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analysis also shows that the relative contribution of the individual 
SSPs, as well as the inferred SFHs (see i]4.2t . are preserved provid- 
ing further confidence on the robustness of our results and method. 



4.2 Inferred Star Formation Histories 

Using full spectral synthesis, we are determining the best-fit combi- 
nation of SSPs at different ages and Zs that make up a given spec- 
trum. As such, the SSPs included and their relative contributions to 
the spectrum can be regarded as a stochastically sampled SFH. We 
can thus examine the radial variation in SFH for each galaxy and 
look for signatures of star formation episodes that have occurred 
at specific locations within the galaxy. Fig. [7] displays a graphical 
representation of the radial variation in SFHs for all eight galax- 
ies in our sample. In each panel, the vertical axis represents galaxy 
radius with the dotted line at r = 0", and the horizontal axis repre- 
sents SSP age. The point size is proportional to the relative weight 
in the fit and the colours and point types code the SSP metallicity. 
Left panels show the light-weighted contributions and right panels 
show their associated mass-weighted contributions. 

Looking at a specific example, examination of the SFH for 
N7490 reveals that the bulge region is dominated in light and mass 
by an old population with solar metallicity. Further out into the 
disk, the light and mass become more dominated by a somewhat 
younger and more metal-poor population. It appears that there was 
an episode of high metallicity SF 1 Gyr ago whose biggest contri- 
bution was at about the bulge/disk transition region (r e = 5.3" for 
N7490), i.e. in a ring at that radius. While the contribution of this 
recent SF episode to the light is quite significant at some radii, it 
does not contribute significantly in terms of mass. Alternatively, 
looking at N7495, which has significant emission in its spectra at 
most radial bins (see Fie. lAlt , we see that the light is dominated 
by a very young (~ 1 Myr), metal-poor population, likely due to a 
recent accretion of pristine gas, but the mass is still dominated at 
all radii by an old and metal-rich population. 

In Table |4] we list for each galaxy the %-light and mass con- 
tributions of "young" (0.001-1 Gyr), "intermediate" (2-7 Gyr), and 
"old" (10-20 Gyr) age SSPs to the fits at three representative radii; 
r = 0, r e , and r e . With only a few exceptions (r<j for N0628 & 
N7741 and r e for 10239), the contribution to the stellar mass by 
the "old" population is > 65% at all radii, while its correspond- 
ing contribution to the light ranges from ~ 10-90%. Again with 
those exceptions, the contribution to the mass by the "young" SPs 
is < 25%, while its light contribution is in the range ~ 10-90%. 

Such observations are important for galaxy evolution studies 
in terms of the quantities that most strongly affect the SP content 
and predictions of a particular model. If one is more interested in 
the dominant component in terms of mass, the current SED can 
be very misleading, especially when interpreted in terms of SSP 
values. As we demonstrate in Appendix IB 3 1 many of these inte- 
grated late-type galaxy spectra are simply not well represented by 
single SSPs, particularly when there is current or recent SF (see 
Fig. IB3I right panel). While the concept of an integrated spectrum 
of a galaxy being dominated in mass by an old populations, but 
in light by a minor "frosting" of recent SF is not new (e.g. Trager 
et al. 2000), the practice of using "SSP-equivalent" ages is still in 
widespread use and could lead to misleading interpretations in their 
comparison with different types of galaxy evolution models. 

Fig.[7]also demonstrates the relatively few SSP templates that 
enter each fit (see Table [3), as noted above. We return to a discus- 
sion of these SFHs in £|5 . II 



Table 4. Percent Contributions of all SSPs in Given Age Ranges 
(Young = 0.001-1 Gyr; Intermediate = 2-7 Gyr; Old = 10-20 Gyr) to Fit for 
Spectra at r = 0, r e , and (when reached) Weighted by Light (V-band 
normalized) and Mass. 
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4.3 Stellar Population Gradients 

Having refined our population synthesis procedure, we can now 
investigate the measured age, metallicity, and dust extinction gra- 
dients in our sample galaxies. Fig. [8] presents the radial profiles 
of the average age (top panels), average Z (middle), and fv (bot- 
tom), for all our galaxies. The solid lines connect the radial bins for 
light-weighted quantities, while the dashed lines are for the corre- 
sponding mass-weighted quantities. The vertical dotted lines indi- 
cate the position of the bulge effective radius and disk scale length 
(when reached). The profiles for both sides of the slit about r = 
are shown, with one side distinguished by solid circles at each bin. 

A first glance of these gradients reveals the remarkable con- 
sistency of the profiles on either side of the galaxy. A perfect match 
would not be expected given that spiral galaxies are not axisym- 
metric systems, but the similarities again provide confidence in the 
stability of our fitting technique for assessing relative trends (as all 
radial bins are modeled independently). 

Any difference between the light- and mass-weighted quanti- 
ties indicates recent SF activity on top of an underlying older SP. 
Indeed, most of the galaxies in our sample show evidence of recent 
SF all the way to the center (the only exception being N7490 which, 
as we will see below, also has the largest central velocity dispersion 
in our sample). Additionally, to the extent that we measure them, 
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Figure 7. Light (left) and mass (right) weighted SFHs for all eight galaxies. In each panel, the vertical axis represents galaxy radius with the dotted line 
at r = 0", and the horizontal axis represents SSP age. Point size is proportional to the relative weight in the fit and the colours and point types code the 
SSP metallicity. Left panels show the light-weighted contributions and right panels show their associated mass-weighted contributions. In terms of mass, the 
integrated spectrum at most radii is generally dominated by old and metal-rich population. The contribution of very recent star formation, while it can dominate 
the light (cf. N7495), is typically only a minor contribution to the stellar mass. 



all disks show evidence for some contribution from an underlying 
old (> 10 Gyr) population (see Table[4]for percentage contributions 
at r d ). 

The dust content is generally small (fy < 3) and is always 
associated with very young SSPs (~ 1 Myr), as would be expected 
since star forming regions are the sites of dust formation. The small 
implied dust extinction may seem at odds with the results of Driver 
et al. (2007) who, based on inclination-dependent bulge and disk 
luminosity functions for a sample of 10095 galaxies, infer face-on 
B-band extinctions of up to 2.6 mag for bulges and 1.1 mag for 
disks. However, if the dust is largely concentrated in the plane of 
the galaxy, as is the case for the models they use with the scale- 
height of the thickest dust component equal to ~ 20% of the bulge 



effective radius, the overall effect for a face-on bulge is that of near- 
complete attenuation of the stellar light behind the dust lane, rather 
than a reddening from light passing through the dust. Thus, even if 
such large attenuations exist in our bulges, because of the face-on 
orientations, we are observing non-extincted light from above the 
dust lane and thus would not expect our spectra to exhibit strong 
signatures of dust reddening. 

The gradients themselves show a wide variety of behaviors. 
Within the bulge-dominated region, the light-weighted age pro- 
files decrease in age with radius for N0173 & N7490, increase 
for N0628 & N7741, and approximately flat for U2124, N7495, 
N7610, & 10239, while the light-weighted Z profiles within the 
bulge decrease with radius for N0173, N0628, & 10239, increase 
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Figure 7 - continued 



for N7495 & N7741, and are approximately flat for U2124, N7490 
& N7610. The disks, on the other hand, almost always show mildly 
decreasing to flat profiles in both age and metallicity. This diversity 
in SP gradients, and the fact that a trend in one is not necessarily 
associated with a trend in the other, allow for a variety of forma- 
tion scenarios. This variety in observed gradients is also in quali- 
tative agreement with those based on Lick-indices of Moorthy & 
Holtzman (2006) for a sample of 38 spirals. We further discuss the 
gradients for individual galaxies in <15.1I also in the context of their 
kinematic profiles which are presented in the following section. 

4.4 Stellar Kinematics 

We now turn to an analysis of the kinematics of our spiral galaxies. 
Given the intricacies involved in accurate measurement of absolute 
velocity dispersion and rotation profiles with long-slit observations, 
we provide below a fairly detailed description of our procedure. 



4.4.1 Measurement of Stellar Velocity Dispersion & Rotation 
Profiles 

For velocity dispersion, a, and radial velocity, V rot , estimates we 
used the Movel algorithm (updated from Gonzalez 1993), first in- 
troduced in <13,4.1I The population templates are now the composite 
stellar population fits obtained in i]4.1l 

In addition to a careful characterization of the resolution and 
wavelength calibration of our GMOS data (§2), in <13.4.1I we also 
explored in detail the absolute spectral resolution, as a function of 
wavelength, as well as the wavelength scale of the BC03 models. 
For the derivation of absolute as and rotation curves, the templates 
were first corrected by the systematic error in their wavelength cali- 
bration derived in <13.4.1l (note that for the stellar population fitting, 
it was the galaxy data that were distorted to the STELIB wavelength 
scale). 

The GMOS resolution is very well characterized by the convo- 
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lution of a top-hat slit resolution of 10.81 ± 0.01 A and a Gaussian 
instrumental resolution of 0.8 ± 0.02 A. Strictly speaking, however, 
this resolution only applies to regions where the galaxy effectively 
fills the 2" slit. Towards the center of each galaxy, where the galaxy 
structures become significantly smaller than the slit width, an ob- 
ject can no longer be considered extended, and the actual resolu- 
tion is given by the convolution of the galaxy surface brightness 
profile, the seeing FWHM, and the spectrograph internal (after the 
slit) resolution (0.8 A). Using the observed spectroscopic bright- 
ness profile of each galaxy and taking into account the measured 
anamorphic amplification, we created 2D effective resolution maps. 
In the Movel procedure, each galaxy is first convolved with the 
residual resolution of the templates (FWHM = \/3.4 2 — 0.8 2 A and 
a = 1 1 km s _1 ) while the template for each galaxy radial bin is con- 
volved by its corresponding row in the 2D-effective-resolution map 
of the galaxy. In this fashion, Movel can then extract absolute ve- 



locity dispersions, since both galaxy and templates are now well 
matched at all points along the slit. 

The templates are generated for bins along the slit that have 
been coadded to accumulate S/N/A ^ 50. To avoid artificially 
broadening each radial bin, the bins must be coadded after remov- 
ing the rotation shape. This is achieved by iteration, where we first 
interpolate the template fits for each coadded bin and use them in 
the unbinned data (continuous 0.07"pixel sampling along the slit) 
to derive and remove the rotation curve. We then reconstruct the 
minimal S/N bins and derive the absolute as and the small residual 
correction of the radial velocity at each coadded bin. 

In Fig. [9] we present the results from our Movel determi- 
nation of the galaxy kinematic profiles. For each galaxy, we plot 
the velocity dispersion profile (top panel), rotation curve (middle 
panel), and 7 profile (bottom panel). Note the different scales for 
each galaxy. Recall that 7 is a measure of the relative absorption 
strengths between galaxy and template, with 7=1 representing a 
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Figure 7 - continued 



perfect match. In addition to obvious potential culprits inherent to 
the data reduction process (sky-line subtraction, relative flux cali- 
bration, etc.), strong deviations from 7=1 could indicate a region 
in parameter space that is not represented in the models. For ex- 
ample, a mismatch of abundance ratios between the models (which 
only include solar-neighborhood abundance patterns) and galaxies 
as described in Sj3AT](see, e.g. U2124), or at low age/Z where the 
stellar libraries are sparse (see, e.g. N7741). It could also indicate a 
spectrum whose emission-line fill-in did not get fully masked out in 
our "cr-clipping" procedure described in ^3.3.11 For the most part, 
however, 7 is quite close to 1 indicating a good match between 
galaxy and model fits. 

The kinematic profiles generally show a rapid decrease in ve- 
locity dispersion outside the bulge effective radius and fairly well 
defined rotation curves. However, given the face-on nature of this 
sample, the correction to absolute rotation speed (i.e. sin(i)) is 
quite large and uncertain. For N0628 and U2124, the odd rotation 



profiles reveal a poor alignment of the slit on the nucleus of the 
galaxy. Regardless, our a profile for N0628 is very well matched 
to that of Ganda et al. (2007; hereafter Ganda07), indicated by the 
red lines, which was derived from 2D SAURON data. 

The dominance of rotational versus random motion support of 
the bulges of our spiral galaxies can be isolated and compared with 
pure elliptical systems in the (V maa; , re /<to, e) plane. Here, V max , rc 
is the maximum rotational velocity measured within the effective 
radius of the bulge, ero is the central velocity dispersion, measured 
within an effective aperture of r e /10 (and labeled as o> e /io, see 
below). The bulge ellipticity, e, also measured at r e , is derived from 
the isophotal SB fitting (see Mac03 for details). The location of our 
bulges in this plane is shown in Fig.[l0] The solid curve represents 
the location of oblate spheroidal systems with isotropic velocity 
dispersions that are flattened only by rotation (Binney & Tremaine 
1987). Three of our bulges lie above the oblate-line, indicative of 
the disky behavior expected of bulges that have undergone secular 
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Figure 9. Kinematic profiles for our eight galaxies. Plotted in each panel are the velocity dispersion, a, in kms" 1 (top), rotational velocity, V ro t, in 
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evolution, whereas a distinction between "pseudo" and "classical" 
for the other four bulges, which lie close to, but below the oblate- 
line, is not as clear-cut. Note that the large uncertainties in N7741 
result largely from the very low velocity dispersion of this galaxy 
and its large ellipticity is due to its dominant stellar bar. 



The main kinematic parameter of interest for this work is the 
central velocity dispersion. This was measured as a light-weighted 
average within a effective aperture of r e /10 for direct comparison 
with the elliptical galaxy sample of Thomas et al. (2005; hereafter 
Thomas05). The average ages and metallicities within this central 
r e /10 aperture were also computed and the values are listed in Ta- 
ble [5] 



4.5 Age and Metallicity vs. Central Velocity Dispersion 



We can now assess any trends in the age and metallicity of our spi- 
ral bulges as a function of velocity dispersion. In Fig. [TJJwe plot 
light- (solid black squares) and mass- (red open squares) weighted 
average age (top panel) and Z (bottom panel) for our eight bulges. 
Simple linear least-squares fits to both trends are shown in the leg- 
ends. The light-weighted values follow the expected trend of in- 
creasing average age and metallicity with larger £>>,,/ 10 ■ However, 
when regarded from a mass-weighted point of view, both trends es- 
sentially disappear, indicating that the mass of all spiral bulges is 
dominated by a population of very old and metal rich stars. 
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Figure 9 - continued 



Table 5. Average Age, Z, and cr within r e /10 



Name ay e/10 do> e/10 (A); (A) m (Z) ; (Z) m 
NGC (kms -1 ) (Gyr) (Gyr) 

(1) (2) (3) (4) (5) (6) (7) 



N0173 


73.01 


5.99 


10.28 


15.47 


0.027 


0.023 


N0628 


55.14 


6.39 


5.33 


11.98 


0.037 


0.028 


U2124 


139.24 


4.98 


10.49 


15.34 


0.029 


0.028 


N7490 


140.92 


5.53 


17.11 


19.19 


0.028 


0.028 


N7495 


96.87 


8.61 


6.57 


14.13 


0.011 


0.032 


N7610 


67.60 


12.57 


7.13 


14.15 


0.016 


0.040 


N7741 


12.17 


13.61 


3.30 


14.41 


0.011 


0.041 


10239 


61.24 


11.80 


4.45 


12.71 


0.041 


0.038 



Notes — Central galaxy parameters averaged within a radius of r e /10 X 2" (the slit-width). Col. 
(1): Galaxy ID. Col. (2): Velocity dispersion within r e /10. Col. (3): Velocity dispersion error within 
r e /10. Col. (4): Average V-band light-weighted age within r e /10. Col. (5): Average V-band 
mass-weighted age within r e /10. Col. (6): Average V-band light-weighted metallicity within 
r e /10. Col. (7): Average V-band mass- weighted metallicity within r e /10. 



4.6 Comparison With Other Studies 

Most studies of the SPs of spiral bulges and ellipticals to date have 
focused on Lick-index measurements. As outlined in the Appendix, 
such measurements are subject to a number of limitations, partic- 
ularly in the case of star-forming systems. The biggest limitation 
is due to emission-line fill-in of the strongest age discriminators, 
the Balmer absorption lines. With high quality spectra, however, 
this limitation can be largely overcome with a careful subtraction 
of the emission-line contamination to the indices, as was done in 
the work of Ganda07 based on SAURON data. The method used is 
described in Sarzi et al. (2006) and was also applied in the study 
of Morelli et al. (2008; hereafter Morelli08). We thus consider the 
age and Z estimates in the Ganda07 and Morelli08 samples to be 
appropriate for direct comparison with the Lick-index-based esti- 
mates for early- type galaxies in the sample of Thomas05. All of 
the above studies provide SSP-equivalent values, and we note that 
these are "light-weighted" and not strictly "average" values in the 
same sense as ours which implicitly measure the fractional contri- 
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bution of each SSP to the galaxy flux, normalized to the V-band. As 
such, small differences in the fractional light contributions of dif- 
ferent SPs can cause large scatter in the derived parameters. This 
effect is clearly seen in Fig. [T2] where we plot the results from our 
study along with those from the early-type sample of Thomas05 
and the bulge samples of Ganda07 and Morelli08. See the figure 
caption for a description of the different data sets. 

A direct comparison of absolute ages is difficult due to differ- 
ences in the models used to derive them. Of particular note are the 
extreme ends of the age ranges used: the models of Thomas et al. 
(2003), that are used in all three of our comparison studies, cover 
SSP of ages 1-15 Gyr, where as the largest age SSP from the BC03 
models is 20 Gyr and we include those as well as 16 Gyr SSPs in 
our modeling (see Fig.[2]l. Thus, our age predictions will tend to bi- 
ased to larger ages whenever an SSP of age > 13 Gyr is included. At 
the young end, we include ages from BC03 down to 1.3 Myr, thus 
our ages could be biased young whenever SSPs of age < 1 Gyr are 
indicated. This would have the overall effect of increasing the slope 
of the age-CTo relation derived from our methods. The black dashed 
line in Fig.[T2]is the best fit scaling relation for low-density envi- 
ronments derived in Thomas05 based on Monte-Carlo simulations 
(which assume linear correlations of the parameters age, metallic- 
ity, and a/Fe ratio with log a and account for observational errors 
and intrinsic scatter in all 3 parameters). Indeed, if we compare 
the slope of the Thomas05 fit with our own derived in Fig. II Hand 
shown again as the black dotted line in Fig. [12] our slope is signif- 
icantly steeper (we note, however, that our statistics are very low). 

Comparison of the age-ao relation (upper panel in Fig. 1 12t. 
for each data set reveals a similar general trend of increasing age 
with (Jo. However, for the early- type sample of Thomas05 (red 
open circles), the scatter in this trend is much larger. This can 
be easily explained in terms of the above-mentioned difference in 
our quoted ages; ours are true light- weighted averages, whereas 
those of Thomas05 are SSP-equivalent values which are heavily 
weighted towards the most recent episode of SF and largely insen- 
sitive to any underlying old population (even though it may domi- 
nate the stellar mass). The same applies for the SSP-equivalent ages 
of the bulge samples of Morelli08 (green open stars) and Ganda07 
(light blue open triangles). However, Ganda07 also attempted to 
model their Lick-indices in the context of an exponential SFH and 
derived the best-fit e-folding timescale, r exp . We converted T exp 
into an average age using Eq. 1 1 in MacArthur et al. (2004) and the 
results are shown in Fig.[l2]as the solid blue triangles. Not only do 
these show much smaller scatter than their SSP-equivalent values, 
but they also follow a very similar trend as our average ages. Thus 
it seems that, when a reasonable attempt to sample the SFH of a 
galaxy is made, the trend of smaller ages with smaller ao is much 
weaker and has significantly less scatter. We speculate that most of 
the scatter in the Thomas05 and Morelli08 samples would disap- 
pear if they were modeled in the context of a realistic SFH (as has 
been previously suggest for pure ellipticals in the "frosting" models 
of, e.g., Serra & Trager 2007). 

All data sets again follow the same, and expected, trend of 
increasing metallicity with ao- The scatter is quite small at high 
velocity dispersion (ao > 100 km s -1 ) but, for the SSP-equivalent 
estimates becomes quite large at smaller ao, again likely due to 
the fact that they are influenced by the stronger recent SF. On the 
other hand, our average values of Z follow the continuation of the 
Thomas05 slope and our derived slopes are almost identical, de- 
spite the small overlap in the ao -range of our respective samples. 

The magenta square in both panels of Fig. [12] represents our 
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to find the best e-folding time, r (see their Fig. 27), where we converted r 
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fit Zs for this method). Green Open Stars: Morelli et al. (2008) cluster SO 
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light-weighted SSP values. Magenta Open Square: Our results for N0628 
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a) is in the right direction in both age and Z, but not quite as strong in Z 
to consider it a true "outlier". The black dashed lines are the best fit scal- 
ing relations for low-density environments derived in Thomas05 based on 
Monte-Carlo simulations (which assume linear correlations of the parame- 
ters age, metallicity, and a/Fe ratio with logcr and account for observational 
errors and intrinsic scatter in all 3 parameters). The black dotted lines are 
the same as in Fig. llll and are least-squares fits to our light-weighted values 
(black squares). Error bars on all parameters for our data are shown on each 
point. The legend in the top left corner of each panel shows the average er- 
rors on the samples from the literature, when provided (Ganda07 do not list 
errors on derived ages and Zs). 



results for N0628 but using the (A); and (Z); measured at ~ 1 r e 
(where the age peaks in our gradients, see Fig.©. The purpose here 
is to compare with the MW "outlier" in Thomas & Davies (2006; 
Fig. 4), as described in [JT] The deviation (from the central value; 
see the black square at the same a) is in the right direction in age, 
i.e. if the age was determined at the bulge effective radius rather 
than at the galaxy center, this galaxy would similarly stand out as 
an outlier in the age - oo relation. We also see evidence of super- 
solar [a/Fe] at these radii in the (Fe) -Mgb Lick-index plane. The 
photometric (Zoccali et al. 2003) and spectroscopic (Rich & Origlia 



2005; Zoccali et al. 2006) studies finding old and a-enhanced SPs 
of the MW bulge indeed sample a larger physical radius, outside 
of the plane, which is different from the "central" values measured 
here for external galaxies. Thus, our observation of a positive age 
gradient, along with the observations of Peletier et al. (2007) of 
multicomponent (classical plus pseudo) bulges, the position of the 
MW can be readily understood and is not a major cause for concern. 
We further discuss the reliability of this result in i]5.1l 

Finally, also of particular note here is that between all three 
samples, a very large range in spheroid mass and type is covered; 
from pure ellipticals, through SO bulges, "classical" spiral bulges, 
and to proto-typical spiral "pseudobulges". However, in the con- 
text of the SP parameters as a function of the central galaxy poten- 
tial, when a proper accounting of the spheroid's SFH is considered, 
there does not seem to be a break in the trends going from pure el- 
lipticals to late-type bulges, nor is there a clear distinction between 
classical and pseudobulges. 



5 DISCUSSION 

Our study is the first to provide radially resolved spectra well into 
galaxy disks, enabling a comparison of the bulge and inner disk 
populations. While the small size of our sample does not allow for 
a detailed statistical study of galaxy parameters with inferred age 

6 Z, we can still infer plausible SF mechanisms with a detailed 
look at each galaxy in the context of specific formation scenarios. 
The currently favored scenarios for galaxy bulge formation can be 
summarized as follows: 

"CLASSICAL" FORMATION 

For the context of the following discussion, we consider "classical" 
formation to imply that the bulges formed on a rapid timescale, as 
is thought to be the case for pure ellipticals. There are two classes 
that fall into this heading as follows; 

Monolithic collapse: galaxies formed early through the 
gravitational collapse of a single cloud of primordial gas (e.g. 
Eggen, Lyden-Bell, & Sandage 1962; Larson 1974; Carlberg 1984; 
Thomas et al. 1999). The relevant predictions from this model in- 
clude: 

• spheroids are metal-rich in center and have a negative metallicity 
gradient 

• spheroids are old and have shallow positive age gradients (centers 
slightly younger) 

• stellar abundance ratios of spheroids are super-solar [a/Fe] with 
positive gradients 

• light profiles of spheroids have large Sersic n values (i.e. high 
concentration) 

• disks accrete onto an already formed spheroid and are thus 
younger 

Merging: galaxies are gradually assembled in a hierarchical 
process through multiple mergers of smaller subgalactic units. This 
process begins at early times and continues to the current epoch 
(e.g. White & Rees 1978; Cole et al. 1994; Thomas et al. 1999; 
De Lucia et al. 2006). 

Predictions for merger scenarios are clearly going to be much 
more complicated and depend on many different factors including 
the specific merger history of a given galaxy, the mass and gas con- 
tent of the infalling clumps at the time of merging, and the specific 
physics describing the SF as a result of the merger and the over- 
all shape of the merger product. In particular, predictions about 
SPs of galaxies, particularly those of multiple-component (bulge, 
disk, bar), are still in their infancy and must currently resort to 
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semi-analytical recipes to describe many of the complicated gastro- 
physical processes that remain poorly understood and too difficult 
to model explicitly. However, certain global predictions based on 
cosmologically-motivated merger trees are robust and include: 

• a galaxy mass-metallicity relationship (more massive galaxies 
have higher-Z) 

• early mergers will have super-solar [a/Fe], but more recent merg- 
ing leads to solar [a/Fe] with flat profiles 

• less massive spheroids have more extended SFHs leading to pos- 
itive correlations of oq with average age, Z, and [a/Fe]. 

In the above two scenarios, the formed spheroid is considered 
a "bulge" if it happens to have formed a disk around it and/or if the 
disk component of a progenitor galaxy survived the merger event. 
Otherwise, the spheroid is considered a pure elliptical galaxy. 

SECULAR EVOLUTION 

In the secular scenario, galaxy bulges formed through a secular re- 
distribution of material from the disk (e.g. Kormendy & Kennicutt 
2004, and references therein). 

The predictions of SP parameters and their gradients in the 
context of secular evolution are also complicated and remain in 
their infancy. A major contributor to the uncertainties involves the 
inclusion of gas in the simulations leading to very different results 
for dissipational versus dissipationless evolution. Furthermore, pro- 
cesses that fall under the guise of secular evolution can occur over a 
wide range of timescales depending on bar presence/strength, spiral 
arm strength, triaxiality of DM halo, etc. 

The more tractable predictions include: 

• dissipationless evolution rearranges disk stars into a bulge com- 
ponent, thus both components would have the same SP parameters 

• dissipational evolution will result in fresh SF in the center leading 
to bulges containing young stars. 



5.1 Notes On Individual Galaxies 

In the context of the above bulge formation scenarios and their pre- 
dictions for bulges properties, we now review each of our eight 
galaxies individually. In Table [6] we provide a summary of the for- 
mation scenarios indicated by all of our observations. 

N0173: The SP populations in the bulge region of this Sc 
galaxy reveal an old and metal-rich center with gradients to 
younger ages and smaller Z within 1 r e , consistent with the "clas- 
sical" monolithic picture of bulge formation. The disk is some- 
what younger with virtually no gradient in age, but a mild gradi- 
ent to smaller Z, out to ~ 2rd- The SFH in Fig.|7]reveals roughly 
equal contributions in light (normalized to the V-band) of an old 
(~ 13Gyr) sub-solar (Z = 0.004) population plus a young (lGyr) 
metal rich (Z = 0.05) population within ~ 1 r e . More recent metal- 
poor SF (possibly from accretion of unenriched gas) is evident in 
the outer region around 1 r e . Throughout, the old population is en- 
tirely dominant in mass (right plot in Fig-EJ- There is evidence for 
slightly enhanced [a/Fe] in the center (Fig. [4}, pointing to short 
SF timescales. This galaxy lies in the "ambiguous" region of the 
{Vmax,re/o~o, e) plane in Fig.[l0] but in general seems consistent 
with a classical monolithic bulge formation scenario. 

N0628: The prototypical grand-design Sc galaxy. Here the 
inner bulge is quite young and metal rich, consistent with recent 
central SF, possibly due to secular funneling of gas to the central 
regions. There is no sign of a bar component in this galaxy, but 
this does not rule out secular processes from, for example, a pre- 
existing bar or the strong spiral structure. The strong spiral dust 
lanes on the concave side of the arms are indicative of shocks from 



the pile up of gas when entering the spiral density wave (e.g. Kor- 
mendy & Kennicutt 2004). The gas loses energy at the shock front 
and sinks to the center resulting in central SF. This process occurs 
on much longer timescales than the action of a bar, but if the spiral 
structure is strong (high pitch angle), such a process will contribute 
to the building of a central component. 

As with N0173, we also see a light-dominant high-Z 1 Gyr 
component on top of an old stellar population. Here, the young 
component has more weight by light and mass in the center, such 
that up to ~ 40% of the mass in the central regions was contributed 
by the SF episode that occurred 1 Gyr ago. At increasing radii, we 
see stronger contribution from much younger and metal-poor SSPs. 
This is indicative of accretion of unenriched material onto the disk. 
However, the main contributor by mass out to ~ lr^ is still from an 
old component. One side of our slit is consistent with a disk age gra- 
dient toward younger ages with increasing radius, consistent with 
inside-out formation. However, there are spikes to old ages and the 
other side of the slit is more consistently at old ages. This is likely 
due to arm-interarm crossings, where the older regions coincide 
with an interarm region where only the underlying oldest SPs are 
present. 

Of particular note is the spike to older ages just beyond ~ 1 r e 
for NGC 628. This is precisely the signature noted in SjJJrequired 
to reconcile the discrepant results for the MW bulge. The fact that 
we observe this signature on both sides of the galaxy (whose spec- 
tra were fit independently) renders confidence that this signature 
is real. Additionally, using SAURON observations, Ganda07 ob- 
served the same trend in age for NGC 628 using a fully independent 
approach. In fact, this trend was first implied by the UV colour gra- 
dients presented in Cornett et al. (1994), from which they inferred 
that the SFH of this galaxy varies significantly as a function of ra- 
dius. Ganda et al. (2006) demonstrated that the stellar kinematics 
are characterized by slow projected rotation, and a mild central dip 
in the velocity dispersion profile, indicating a cold central region, 
which could be identified as an inner disk. Our velocity dispersion 
profile agrees very well with theirs, as shown by the red lines in the 
top panel of the top right plot in Fig. [9] 

As mentioned above, images reveal clear spiral dust lanes 
throughout the disk of this galaxy. From the bottom panel in the 
top right plot in Fig. [8] we see that indeed, our fits indicate the pres- 
ence of modest amounts of dust at various radii. 

In all, there seems to be a significant contribution from secular 
processes to the growth of this bulge. 

U2124: This is our earliest-type, strongly barred, spiral, 
SB(r)a; our slit was aligned along the bar. Both the light- and mass- 
weighted age profiles are very similar; quite flat at old ages and 
super-solar Z within 1 r e . The flat profiles could be interpreted in 
the context of redistribution of material due to the bar. However, 
the small weight of the younger SPs indicates that perhaps the gas 
supply is feeble such that only a small amount of SF is triggered by 
the bar action. There is no evidence for dust extinction in the spec- 
tra, which is also consistent with no recent SF and little gas con- 
tent. There is also a clear gradient in Z beyond 1 r e which would 
be consistent with a bulge formation scenario similar to ellipticals 
(i.e. early and short timescale of SF, also indicated by the super- 
solar [a/Fe] in Fig. |4). Unfortunately, our data do not extend past 
the bar radius (~ 21"), so we cannot examine the disk profile. 

The SFH indicates a dominant old population that is a mix- 
ture of a 13 Gyr Z = 0.004 SSP and a ~ 16 Gyr solar/super-solar 
SSP. On top of this, within ~ 1 r e is a metal-rich 1 Gyr population 
extending all the way into the center. This young population con- 
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tributes ~ 50% in V-band light-weight, but only of order 5% to the 
total stellar mass. 

In all, this bulge shows many classical formation features 
alongside pseudobulge-type structures, thus likely is composed of 
a mixture of both bulge types. 

N7490: This is a normal Sbc galaxy. The spectra are quite red 
in the center with weak emission seen in [N II], and a general bluing 
to larger radii. The age profile reveals a very old ~solar Z center 
with a strong gradient to younger ages out to ~ 1 r e and a weaker 
gradient in the disk to young ages out to ~ 2 rd, but with a flat-to- 
increasing Z within 1 r e and decreasing in the disk. These inner SP 
gradients are somewhat ambiguous in terms of their implication for 
the formation of the bulge, but the profiles are generally consistent 
with a merger scenario for the bulge with subsequent inside-out 
disk formation. 

The SFH reveals that the central ~ 1/2 r e is entirely dominated 
(in mass and light) by an old solar metallicity SSP, and this is the 
only galaxy in our sample that has almost no contribution from a 
younger SP in the very central region (to ~ 1/2 r e ). It also has the 
highest (Jo of the sample. This would also be consistent with forma- 
tion by mergers at early times and/or the merger of gas-poor sub- 
components for the bulge component. This is further indicated by 
the super-solar abundance ratio in the center in Fig.|4]as well as the 
strong features in the fit residuals in Fig. [3] similar to those shown 
in Fig. [5] Beyond ~l/2r e , there is a slightly younger (~ lOGyr) 
and more metal-poor population which is likely the accretion of 
a disk after the initial bulge collapse. The SFH also reveals a SF 
episode 1 Gyr ago in a ring with a radius of about 1 r e , where we 
also see the ages dip to lower values and Z to slightly higher. Such 
an episode of recent, enhanced, SF in a ring at 1 r e is a secular pro- 
cess. We also see from Fig. [8] that there is little to no dust in the 
central regions, but significant a "spike" at ~ 1 r<j associated with 
the very young SSP contributing significantly to the light in this 
region. 

This bulge is most consistent with having formed "classically" 
via mergers at early times with subsequent disk accretion, with a 
minor, more recent, secular contribution in the form of a ring of SF. 

N7495: This SABc galaxy has strong emission throughout, 
being strongest at the very center (Fig. lAlb . The current SF dom- 
inates the light, but is entirely insignificant by mass (~ 1-2%), 
leaving an old and metal-rich population dominant by mass out to 
~ 1 rd- There is a significant amount of dust out to ~ 1 r e , again 
associated with the current SF. 

The age profile in the bulge is a bit erratic, jumping from old 
to young and back to old within 1 r e , while Z shows a consistent 
rise in this region. The disk profiles are both very flat and have 
significant contributions to the light of a very young & metal-poor 
SP. Significant amounts of dust throughout are evident, again in 
conjunction with the fresh SF. In mass-weight, however, the disk 
is consistently very old and metal-rich. This is difficult to interpret 
in terms of an inside-out scenario, but given the current level of SF 
throughout, this disk is still in formation. 

There is some evidence of recurrent episodes of SF in the very 
central regions; one contributing -~30% by mass 10 Gyr ago and a 
less massive event (~1%) at 0.4 Gyr; given the stochastic sampling 
of our SFHs, this could be interpreted as a more extended SFH. 

The kinematics reveal a sharp decline in the velocity disper- 
sion from ~ lOOkms -1 at the center, to < 20 kms -1 by 1 r e . Its 
position in the V/ffo diagram (Fig.llOb reveals a significant amount 
of support against collapse from rotation, implying a disky origin, 
whereas the bulge profile shape has our highest value of the Sersic 
n shape parameter, pointing to a more classical bulge. 



Again, there seems to be contributions from both classical and 
secular components to the formation history of the bulge of this 
galaxy. 

N7610: This galaxy is our latest-type class, Scd. There is sig- 
nificant SF in the center of this galaxy. The bulge light profile is 
characterized by a Sersic n = 0.8, highly suggestive of a disky ori- 
gin. The SP profiles are flat within the bulge region with a mild 
trend to younger & more metal-rich inward (which is consistent 
with rapid E-like formation). Again, the old population represents 
the bulk of the mass, but there is significant recent sub-solar SF 
(ages 0.001-1 Gyr). Our data do not extend far enough to asses the 
disk profile, but they appear relatively flat out to several bulge radii. 
Small levels of dust are seen coincident with the regions of strong 
current SF. 

The kinematics reveal a sharp decline in the velocity disper- 
sion from ~70kms~ 1 at the center, to < 10 kms" 1 by 1 r e . Its 
position in the V/ oq diagram in Fig. [To] is below the oblate line, 
but not far enough to place it unambiguously as occupying the E- 
type location of the plot. 

There is evidence for contributions from both classical and 
secular components to the formation history of the bulge of this 
galaxy. 

N7741: This is our only other strongly barred galaxy and our 
latest-type, SB(s)cd. The light-weighted bulge is very young with 
sub-solar metallicity, but in mass-weight is old & metal-rich. The 
recent central SF indicated by the light-weighted values could be a 
result of in-funneling of fresh gas due to the bar. Our slit is aligned 
perpendicular to the bar, thus no flattening of gradients would nec- 
essarily be expected. The disk profile is generally young and solar 
Z with close to flat profiles. 

While there are clear secular processes contributing to the cen- 
tral growth of this galaxy, the overwhelming dominance by mass is 
that of a very old and metal rich whose gradients are consistent with 
a classical picture. 

10239: An SAB(rs)cd galaxy that also shows a significant 
population of 1 Gyr metal-rich stars on top of an old metal-rich 
population. In this case both contribute significantly to the stellar 
mass. The age gradient in the bulge region is relatively flat out to 
~ 1 rd, while the Z profile is decreasing. Disk gradients are gen- 
erally flat to mildly decreasing. There is an indication of slightly 
super-solar [a/Fe] (Fig.0, indicative of rapid formation, while the 
Sersic n = 0.9 light profile suggests a disky origin. 

Once again, this bulge displays a mix of classical and secular 
evolutionary processes. 



6 SUMMARY 

Using radially resolved long-slit spectra of eight star-forming spiral 
galaxies, we have performed a detailed analysis of their stellar pop- 
ulation and kinematic profiles and provide interpretations based on 
currently favored formation scenarios. Central correlations in ages 
and metallicities with velocity dispersion are compared to other 
studies of bulges and pure elliptical galaxies. 

The most pertinent observations are as follows: 

Fitting Techniques: 

• With moderate spectral resolution, good A coverage, and 
high S/N/A (^50), measurement of light- weighted average ages & 
metallicities for star-forming galaxies is feasible. 

• Details are critical: calibration (A & relative flux), resolution, 



© 2008 RAS, MNRAS Q00.[Tlt38l 



Population Synthesis of Spiral Galaxies 25 



Table 6. Summary of "Secular/Pseudo" (S) vs. "Classical" (C) Indicators for Each Galaxy 



Name 


Sersic n 


Bar 


Young SP 


Old SP 


d(A)/dr 


d{Z)/dr 


V/a 


a/Fe 


Overall Type 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


N0173 


S/C 




S 


C 


C 


c 


C 


C 


c 


N0628 


s 




S 


c 


S/C 


S/C 


s 


S 


s 


U2124 


S/C 


S 


S 


c 


S/C 


S/C 


c 


c 


c/s 


N7490 


S/C 




S 


c 


c/s 


S/C 


s/c 


S/C 


c/s 


N7495 


C 




S 


c 


S/C 


S/C 


s 


s 


S/C 


N7610 


s 




S 


c 


c/s 


c/s 


c 


s 


s/c 


N7741 


s 


S 


S 


c 


S/C 


S/C 


s 


s 


s 


10239 


s 




S 


c 


S/C 


S/C 


S/C 


S/C 


s/c 



/Vofes — Col. (1): Galaxy ID. Col. (2): The dividing line in at n~2. Col. (3): If a bar is clearly present. 
Col. (4): Significant presence of recent SF (within the last ~ 1-2 Gyr). Col. (5): Significant presence of very 
old SP (> 13 Gyr). Col. (6): Indication from SP age gradient. Col. (7): Indication from SP Z gradient. Col. 
(8): Indication from kinematics (support from rotation vs. random motions). Col. (9): SF timescale implied 
by central abundance ratio (as inferred from the results of Fig. |4), whereby enhanced oTFe implies a short 
formation timescale (as the Fe elements produced by Type I supernovae have not had time to enrich the gas in 
Fe before cessation of SF). Col. (10): Overall impression. 



velocity dispersion, and rotation must be treated self-consistently 
between the data and models. 

• Different fitting techniques weigh age, metallicity, and abun- 
dance ratios differently: Balmer emission limits age fitting from 
indices; age information is recovered (in the presence of emission) 
from full spectrum and continuum SED fitting (but compounds the 
caveats about dust extinction and fluxing that are less important for 
indices). 

• SSPs are not a good match to late-type galaxies. The degen- 
eracies between age, metallicity, dust, etc., are extreme leading to 
unstable fits. 

• Full population synthesis is the only method to provide reli- 
able and consistent results, and is thus the method of choice for 
the establishment of light and/or mass-weighted average ages and 
metallicities of late-type galaxies. 

Stellar Populations of Bulges and Inner Disks: 

• In a mass-weighted context, all bulges are predominantly 
composed of old and metal rich SPs (> 80% by mass). 

• Some contribution to bulge growth by secular evolution is 
clearly evident in most late-type bulges, with SPs of 0.001-1 Gyr 
contributing as much as 70% to the optical light. The corresponding 
contribution to the bulge mass, however, is generally small (~ 20% 
by mass or less). 

• Spiral bulges display a wide variety of age and metallicity gra- 
dients (from negative to positive) in the bulge region, allowing for 
a range in formation mechanisms. 

• The observation of positive age gradients within the effective 
radius of some late-type bulges helps reconcile the long-standing 
discrepancy of the secular-like kinematics and light profile shape 
(including the presence of a bar) with "classical"-like old and a- 
element enhanced SPs observed in the Milky Way bulge as being 
due to SP sampling at different physical locations in the bulges. 

• Spiral disks show mildly decreasing to flat profiles in both age 
and metallicity, generally consistent with inside-out formation. 

• Bulges follow the correlations of increasing light-weighted 
age and Z with central velocity dispersion as for elliptical galax- 
ies and early-type bulges found in other studies, but when a SFH 
more complex and realistic than an SSP-equivalent is taken into 
account, the trend is shallower and its scatter is much reduced. 

The implication seems to be that bulge formation has been 



dominated by processes that are common to all spheroids, whether 
or not they currently reside in disks. Monolithic collapse cannot be 
ruled out in all cases, but merging must be invoked for most of our 
spiral bulges. The process of formation occurs on shorter timescales 
for spheroids with the highest central velocity dispersions, and the 
relative contribution to the stellar mass budget in bulges via secular 
processes or "rejuvenated" star formation is small, but generally 
increases in weight with decreasing ao- 

Ultimately, we desire a large enough sample to solidify the 
conclusions inferred in this analysis and to assess any trends with 
galaxy parameters, but already these results provide important 
clues for bulge formation scenarios and restrictions for future im- 
plementations of galaxy formation models. 
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APPENDIX A: ID SPECTRA 

In Fig. lAll we present the fully calibrated and velocity (redshift + 
rotation) subtracted spectra at a number of radii for each galaxy. 
The spectra are plotted on a logarithmic scale and adjusted with an 
arbitrary constant to offset the spectra from each other. Only one 
side of the slit is shown, and not all radial bins are plotted for the 
well sampled galaxies, but the r = and last two bins are always 
plotted. The spectra were coadded (radially) to a minimum S/N/A 
^ 50 measured in the 5050-5450 A wavelength interval (selected 
to avoid prominent emission lines such as [O III] AA5007,4959, 
Ha A6563, [N II] AA6583,6548, and the [O I] sky line at A5577). 
The minimum S/N/A of 50 was selected as a compromise between 
accurate age and metallicity determinations, and sufficient radial 
extent for gradient measurements. The light- weighted radius and 
the spectrum's corresponding S/N/A are labeled at the right edge of 
each panel (note that the r = 0" bin is in reality closer to r = 0704 
due to our binning of 072/pixel and assuming a Gaussian light pro- 
file in the galaxy center). Many of the Lick indices (see £]B2t are 
marked as vertical shaded lines - dark gray delineates the central 
passband and light gray the pseudo-continua. 

A number of pertinent observations can be made from a visual 
examination of the spectra: 

• Improper subtraction of prominent sky emission lines (falling 
in and around the NaD and TiO indices) is evident at the largest 
radial bins for all galaxies. 

• There is a significant overall bluing of the SED with radius 
for most galaxies. The two notable exceptions are N7741 & 10239. 
The former shows significant emission in Ha, H/3, & H7A (the 
latter two appearing as fill-in to the underlying absorption line), 
and at [N II] AA6583,6548 and [O III] A5007 in the central spec- 
trum, which could mask an underlying redder stellar SED. N7495 
& N7610 also show significant central emission and bluer SEDs 
compared to those galaxies with less obvious central emission (e.g. 
N0628 & U2124), but a significant (further) radial bluing is still 
evident. 

• For some galaxies with significant emission in the Ha+[N II] 
region, an emission spike superimposed on an underlying absorp- 
tion can be seen in the other Balmer lines (most noticeably in H/3, 
but a spike can be seen even in H7 & H<5 in extreme cases, e.g. 
N7495 & N7741). Given that the Balmer lines are the most sensi- 
tive age indicators in the Lick index system, it is likely that index- 
based age estimates of the underlying stellar population will be dif- 
ficult to constrain in these galaxies. In the full-spectrum fits, these 
emission lines will need to be masked out as emission lines are not 
included in the BC03 models. 

The most likely source of the emission is due to normal, but in- 
tense, current star formation. We can rule out significant contribu- 



tions from active galactic nuclei (AGN) as the source of emission 
based on the low [N II] A6583/Ha and [O III] A5007/H,3 emission 
line ratios (see, e.g., Hao et al. 2005), but small contributions in 
our most severe cases, N7495, N7610, & especially N7741, cannot 
be entirely ruled out (but these would be restricted to Seyfert 2 or 
LINER AGN as the broad lines characteristic of Seyfert 1 spectra 
are not seen). 

• Clear radial variations in many Balmer and metal-line absorp- 
tion line strengths can be seen in many of our galaxies. The inter- 
pretation of these gradients in terms of age and metallicity require 
an implementation of one of the fitting techniques described below 
in §B\ 



APPENDIX B: AGE AND METALLICITY FITS 

Several methods exist to assign an age and metallicity to an ob- 
served integrated galaxy spectrum. Each technique has its own ad- 
vantages and drawbacks depending on the application. One of the 
fundamental challenges any successful method must overcome is 
the inherent degeneracy in stellar populations between old age, high 
metallicity, and dust content, all leading to redder SEDs. Further 
difficulties can come from non-stellar sources that are not included 
in the SP models, such as emission lines from H II regions, and 
contributions from a central AGN. Below we consider 3 methods 
for deriving the SP parameters from our long-slit spectra of star- 
forming spiral galaxies: (i) the Lick/IDS-index system; (ii) SSP fits 
to the entire spectrum; and (iii) full synthesis fits to the entire spec- 
trum. The latter two also include a prescription for dust extinction. 
In the following, after discussing the importance of understanding 
the resolution of both the data and models when making compar- 
isons, we describe the first two methods individually and then com- 
pare their results with each other and with those obtained from the 
full population synthesis fits in §3. 

Bl Resolution and Velocity Dispersion Effects 

In order to compare the data to the models, potential resolution 
effects must be accounted for. In the case where the models have 
higher resolution than the data, this can be accomplished by degrad- 
ing the models to the resolution of the observations. Alternatively, 
when measuring absorption-line pseudo-equivalent widths, the in- 
dices could be corrected for resolution effects using a calibration 
based on a systematic broadening of templates of similar spectral 
type as the data (e.g. PS02). The BC03 models used for this anal- 
ysis have a reported resolution FWHMgo„ ss = 3 A that is constant 
with A. However, as mentioned above in ^3 .41 due to an error in 
accounting for the relative velocities of the stars before combining 
them into SSPs, the effective resolution is not only a bit higher, but 
is best characterized by a constant FWHM with A term plus a ve- 
locity (constant with A A/ A) term. The precise characterization will 
also be slightly different for each model SSP, depending on the rel- 
ative velocity shifts and weights of the stars entering each SSP. On 
average, we find that the FWHM is closer to 4 A, so it is still much 
smaller than the 10.8 A resolution of our data, which is dominated 
by the slit width. Thus, to account for resolution effects due to the 
slit, we convolve all SSP models with a boxcar of width 10.8 A. 

Another broadening effect, independent of the instrumental 
setup, that must be considered prior to any direct model/data com- 
parison is the galaxy velocity dispersion, a v d, along the line-of- 
sight. Galaxy disks are dynamically cold systems, flattened and 
supported by rotation, and thus have small velocity dispersions 
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Figure Al. Spectra as a function of radius for all galaxies. The spectra were coadded (radially) to a minimum S/N/A of 50. Not all radial bins are shown, but 
the first bin and last two bins are always plotted. The light-weighted radius and the spectrum's corresponding S/N are labeled at the edge of the plot. Many of 
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(of order ~ 10-20 km s -1 ), whereas the spheroidal bulges of spiral 
galaxies can have significant support from random motions. Proctor 
& Sansom (2002) measured velocity dispersions for late-type spi- 
ral bulges in the range 50-200 kms -1 for Hubble types Sa-Sbc. 
Because our galaxies are mostly of later type (Sbc - Scd, with one 
Sa bulge), we can expect bulge velocity dispersions closer to the 
small end of this range. 

We have measured the velocity dispersion at each radial bin 
(see i]4.4t for our sample and, in principle, a unique set of mod- 
els could be created for direct comparison for each radial bin of 
each galaxy. However, given that we are analyzing thousands of in- 
dividual spectra, constructing a model for each one becomes too 
onerous. As an alternative solution, we have opted to degrade the 
resolution of all galaxy spectra to that with the highest dispersion, 
so one set of models, properly convolved to this same velocity, can 
be used for direct comparison with the data. We thus account for 
resolution by putting both the data and models on the same scale 
by convolving all observed spectra to a = 168 km s~\ and all mod- 
els with a boxcar of width 10.8 A (corresponding to the slit width) 
and to a= 168 kms -1 . 

B2 Age and Metallicity Fits from Lick Indices 

The Lick/IDS system of spectral line indices was designed to cali- 
brate the strength of fundamental spectral features in stars and com- 
posite systems (e.g., Gorgas et al. 1993). The indices measure the 
strength of a particular spectral feature (either atomic and denned 
as an equivalent width in A, or molecular and measured in magni- 
tudes) relative to a pseudo-continuum on each side of the feature. 
The most reliable indices have been calibrated as a function of stel- 
lar colour (effective temperature), surface gravity, and metallicity 
(Gorgas et al. 1993; Worthey et al. 1994) allowing for the construc- 
tion of semi-empirical population models (e.g. Worthey 1994). 
Each Lick index is sensitive to the metallicity and age of stellar 
populations to varying degrees. When compared with population 
models, diagnostic plots of age versus metallicity-sensitive indices, 
such as H/3 versus Mgb or Fe, help break the age-metallicity de- 
generacy. However, measurements of many of the Lick indices are 
quite sensitive to spectral resolution and, thus, to the velocity dis- 
persion of the system (Gonzalez 1993; Trager et al. 1998; Proctor & 
Sansom 2002), and their use requires relatively high signal-to-noise 
data (S/N/A ^ 50; see Cardiel et al. 1998). As mentioned above, we 
account for the former by convolving all data and model spectra to 
the same resolution before measuring the indices, and the latter by 
coadding the radial spectral bins to S/N/A 50. In addition, certain 
indices, H/3 in particular, can suffer from nebular emission contam- 
ination, even in early-type galaxies (Gonzalez 1993; de Zeeuw et al. 
2002; Caldwell et al. 2003). 

To help overcome the problem of nebular emission fill-in of 
the H/3 feature, Worthey & Ottaviani (1997, hereafter W097) in- 
troduced two pairs of indices that measure the higher order Balmer 
lines H7 and H<5. While their age-sensitivity is not as strong as 
for H/3, the higher order Balmer lines are much less affected by 
emission from ionized gases (e.g. Osterbrock 1989; Osterbrock & 
Ferland 2006). Thus, when combined with a metallicity-sensitive 
index, the W097 indices provide a more reliable age estimate for 
star-forming galaxies. 

Finally, another more subtle issue is the sensitivity of many 
of the indices to non-solar(neighborhood) abundance ratios (e.g. 
TMB03). The BC03 models used here do not account for variations 
in [a/Fe] which are expected for SFHs with different timescales. 
Higher [a/Fe] is associated with very short SF timescales as there 



is not enough time for the Fe produced in Type la supemovae to 
enrich the star forming gas that is heavy with a elements produced 
in the short-lived massive stars that end their lives as Type II super- 
novae. On the other hand, for more extended SF, there is plenty of 
time for Fe to enrich the gas to solar levels of [a/Fe]. We attempt to 
address all of these issues in the context of our spiral galaxy spectra 
in the following. 

In all, there are currently 25 Lick/IDS indices defined and cali- 
brate^. Our spectra cover 24 of the indices (missing only the bluest 
high-order Balmer HSa index). We measure all 24 and use various 
combinations thereof in an attempt to assign reliable light-weighted 
age and metallicities to each spectrum. 

The obvious advantage of using the Lick/IDS-index system is 
that the narrow wavelength span of each index measurement ren- 
ders them less sensitive to low-frequency effects such as overall 
flux calibration errors and dust reddening (see MacArthur 2005). 
One major drawback when considering galaxies with current or re- 
cent SF is that the system was inherently designed to study the old 
stellar populations associated with elliptical galaxies and globular 
clusters. 

The age and metallicity fits based on Lick indices use same 
maximum likelihood approach as in MacArthur (2005), but here 
the method has been extended to include more than just two in- 
dices. So, in the figure of merit (Eq. [A5] in MacArthur 2005), N 
is now an arbitrary number of indices (at least 2 and up to the 24 
we measure here). 

While the errors on the measured indices are independent, the 
model tracks are not orthogonal. This produces non-orthogonal er- 
rors in age and metallicity, which can lead to spurious correlations 
if not understood. Additionally, the degree of non-orthogonality of 
the age and metallicity tracks changes with position on the grid. For 
this reason we have used Monte Carlo methods to model the "effec- 
tive ranges" of the fitted ages and metallicities, taking into account 
a normal distribution about the measured errors (see MacArthur 
et al. 2004 for a more detailed description of this approach). 

We have seen in [jA]that many of our galaxies have signifi- 
cant amounts of emission, most likely coming from star forming 
H II regions. Many authors attempt to correct for emission by fit- 
ting emission-free templates to their galaxy spectra. The templates 
can either be from linear combinations of stellar templates (e.g. 
Gonzalez 1993), or from models (such as the BC03 models used 
here). Neither method is ideal. The former requires a library of stel- 
lar templates, ideally taken in the same observing run with identical 
conditions, that match the galaxian spectra extremely well. The lat- 
ter imposes a model dependence since, if a model with a given age 
is used to make the correction, the same model age will be returned 
(you get out what you put in). Additionally, a significant amount 
of dust in the galaxy would cause a reddening of the spectrum that 
must be included as another model-dependent free parameter in the 
template fits (for both the empirical and model methods). 

In a study of 40 elliptical galaxies, Gonzalez (1993) cre- 
ated individual templates from stellar spectra obtained with the 
same observational set-up for each galaxy spectrum (at each ra- 
dius) for the purpose of measuring accurate velocity dispersions. 
These templates fit the observed spectra to within ~ 1%. Divi- 
sion of the galaxy spectrum by the best fit template enabled de- 
tection of faint levels of emission (seen predominantly in [O III] 
and H/3) in over 60% of his sample. From the galaxy/template 



7 For up-to-date definitions of all 25 Lick/IDS indices, see Guy Worthey's 
webpage at http://astro.wsu.edu/worthey/html/index.table.html 
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spectrum, he computed pseudo-EWs for the [O III] A 5007 & 4959 
A and H/3 emission and found them to be strongly correlated as 
EW(H/3) = 0.7EW([O III]). As such, the H/3 index in absorption 
could be corrected for emission line fill-in by adding the correction 
0.7EW([O III]) to the measured index. He also made a correction 
to the Fe5015 index as the [O III] A5007 line lies within the limits 
of the central bandpass. Gonzalez (1993) strongly emphasized the 
importance of having suitable templates for reliable measurements 
of both velocity dispersion and emission corrections. 

The tight correlation between [O III] and H/3 emission was 
tested for later-type bulges (SO-Sb bulges) by PS02 using a 
method similar to that of Gonzalez (1993). While they found the 
same [O III] -H/3 relation for their earliest-type bulges, the six 
late-type bulges did not follow the Gonzalez (1993) correlation for 
ellipticals, with some bulges showing a stronger H/3/[0 III] ratio, 
while others showed clear [O III] emission, but none in H/3. This is 
likely due in part to the stronger absorption in the younger stellar 
populations, but could also be due to template mismatch, or true 
physical differences in the emission line regions for these galaxies. 
Thus, no straight-forward correction seems to apply to spiral bulges 
and, by extension, their disks as well. 

The higher-order Balmer lines are much less affected by emis- 
sion line fill-in. For a broad range of physical conditions expected 
in galactic emission-line regions, the line-strength ratios of H5 & 
H7 relative to H/3 are 0.25 and 0.5, respectively (Osterbrock 1989), 
providing some relief from emission line fill-in. Indeed, in many 
of our galaxies, the H/3 index lies far off of the model grids (from 
which ages ^> 20 Gyr would be inferred), whereas the higher-order 
indices imply much younger ages (~ 5 Gyr). However, it is still 
possible that the higher-order Balmer lines may suffer significant 
fill-in in strongly star forming regions (and indeed we see this in 
a few of our spectra; see, e.g. the central spectrum of N7495 in 
Fig.HD. 

Given the above mentioned issues about emission-line correc- 
tions to the measured Balmer-line indices, we have instead opted 
for scheme whereby the emission line affected indices are system- 
atically eliminated from the fit, and the age determinations are com- 
pared. We prefer this method as it is less model-dependent and 
takes full advantage of the narrow baselines of the indices, while 
still using information from the entire observed SED. 

We also attempt to address the issue of abundance ratio vari- 
ations. While predictions for the effects of enhanced [a/Fe] exist 
for Lick indices (e.g. Trager et al. 2000; TMB03; Thomas, Maras- 
ton, & Korn 2004; Schiavon 2007), these do not yet provide high- 
resolution SSP spectra, thus we cannot properly account for res- 
olution by matching the models to the data, nor can we do full 
spectrum fits for comparison. To date, the only models that provide 
high-resolution spectra for stars with non-solar abundance ratios 
are those of Coelho et al. (2007). Comparing the results obtained 
here with these new models will be the focus of an upcoming paper, 
but for now, we consider only the BC03 solar-scaled models for all 
fits. 

Our fitting scheme for deriving light-weighted SSP ages and 
metallicities from Lick indices considers 5 different cases as fol- 
lows: 

Case A: Fit All Well-Measured Indices 

By "well-measured" we mean that the index (all 3 bands) is 
not severely affected by systematic problems such as the wave- 
length gaps between the 3 GMOS CCDs or strong sky lines (which 
are extremely difficult to subtract accurately). The indices that were 



systematically eliminated from all fits are highlighted in gray in the 
leftmost column in the deviation plots (left panels in Figs. lBlHB2l 
see below for plot description). Note that we never include the NaD 
index in any of the fits as this index is affected by sky lines as well 
as poorly understood absorption from the interstellar medium of 
the galaxy. The TiO indices were also excluded from all fits except 
for N0628 due to sky line subtraction issues. 

Case B: Case A (All Balmer-Line Indices, Fe5015, & CNi) 

Here we attempt to measure ages and metallicities by eliminating 
all indices severely affected by emission. This is feasible because 
of the significant age dependence of G4300 and the weak age de- 
pendence of the metal-line indices. It is clear that all Balmer-lines 
could suffer from emission line fill-in. The Fe5015 index has the 
[O III] A5007 line in its central bandpass thus could also suffer fill- 
in. While seldom mentioned in the literature, we note that the CNi 
index could also be compromised by a strong source of emission 
since it contains the H<5 line in its blue continuum (which forces the 
index to larger values). This is clearly the case for at least two of 
our galaxies, N7495 & N7741, where H<5 is seen in emission (see 
their respective spectra in Fig. IAU . As an example of the strong 
influence that emission can have on CNi, including CNi in the 
fits, having removed all Balmer-line indices and Fe5015, the age is 
forced to the maximum model age of 20 Gyr to accommodate the 
large CNi . When CNi is removed, the best-fit age drops to 1-2 Gyr. 
We also remark that, because of our resolution, if there is a strong 
emission spike in H7, its tails can bleed into the red passband of 
the G4300 index, but G4300 was not eliminated from the fits. 

Case C: Case A - (H/3 & Fe5015) 

Eliminate only the indices most severely affected by emission. This 
case should be appropriate for galaxies with only small amounts of 
emission, which could be undetectable in the spectra, because the 
higher-order Balmer-lines will be only weakly affected. If the age 
estimates here are younger than those in Case A, emission is likely 
present, whether or not it is clear from the spectra. 

Case D: Fit Only Hy A Hy F , [MgFe]', & G4300 

These indices are the least affected by abundance ratios (TMB03; 
Thomas, Maraston, & Kom 2004). While only a few of our galax- 
ies showed evidence for significant [a/Fe] enhancement (N0173, 
U2124, & N7490; see Fig.0}, many showed CN 2 values far be- 
yond of the model grids (N0628, N7490, 10239), which this com- 
bination of indices also avoids. Such CN enhancements have been 
observed in globular clusters and elliptical galaxies (Henry & 
Worthey 1999), and are likely due to an enhancement in nitro- 
gen (rather than carbon). However, the interpretation of carbon- 
enhancement is challenging due to the many, and poorly under- 
stood, production sites for carbon. 

Case Z: a Clipping of High Deviators 

Here we initiate the fit using all "well-measured" indices. We 
then follow an iterative procedure removing at each step indices 
that have large data— fit deviations: \xi\ = |data; — fiti|/<5i > Xmax, 
where Xi is the relative data/model fit deviation for index i in units 
of the error on the measured index, 8%. The threshold Xmax is set 
to 3.5 initially, and reduced by 0.5 with each iteration until a mini- 
mum value of 1. Thus, all indices with Xi ^ 1 get eliminated from 
the fit, and ^ 2 indices must remain for a successful fit. 
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Results from the age and Z fits for all five cases are displayed 
in Figs. lBTHB2l for U2124 and N7495, respectively. In the left pan- 
els of each plot, we show the deviation in units of error for each 
index in all four fits, quantified as \i = (data* — fiti)/<5,. Each case 
is denoted by a different colour: Case A: green, Case B: blue, Case 
C: pink, Case D: orange, and Case Z: red. Three-pronged point- 
types mark the indices included in the fit for each case, while open 
squares mark indices omitted from the fit. Point sizes and rotations 
were adjusted for visibility when points are overlapping. Thus, if a 
given colour is not visible for a given index (except (Fe), which is 
never explicitly fit, and [MgFe]' for which we only show Case D 
results), it is because the \ for that index is off scale. The dotted 
black vertical line mark zero deviation. The indices are labeled on 
the y-axis (the dashed horizontal lines are there to guide the eye 
to the index labels); those marked in gray are the poorly-measured 
indices that are never included in the fits (note they are different for 
each galaxy). The index labels on the left indicate those excluded 
from all fits, while the labels on the right indicate the indices that 
ended up included (blue) in the Case Z fit for the r = 0" spectrum. 
For each galaxy we plot deviations at three different radii (labeled 
in green above each panel): the central point [left panel], the point 
closest to the bulge r e [middle panel], and the point closest to the 
disk scale length r& [right panel]. Below the radius label for each 
panel are the best fit ages and metallicities for each case in their 
associated colours. 

The right panels in Figs, IBlHB2l show the age gradients in log 
[top panel], metallicity gradients in log [middle panel] and the \ 2 
figure of merit [bottom panel ], as a function of the logarithmic ra- 
dius scaled to the disk scale length (with labels in arcsec on the top 
panel) for all five fit cases (the colour scheme is the same as in the 
deviation plots and is labeled at the bottom of each plot). The data 
points for each fit are connected by dotted lines to guide the eye. 
Labels on the right axes show the linear values. For each galaxy, the 
seeing FWHM, s, the bulge r e , and the disk are indicated with 
arrows at the top. The dashed lines in the age and metallicity plots 
denote the model limits, so fits with these values should be treated 
with caution (recall that we do not extrapolate beyond the model 
grids). 

Examining the fits for the different cases in Figs. lBll4B2l we 
see that the different cases can yield very different best-fit parame- 
ters. The reason for this in some cases is clear, for example, when 
there is significant emission and the fits are different depending on 
which of the Balmer-lines (if any) are included in the fit (see, e.g., 
N7495 fits in Fig. lB2t . From the deviation plots we also see that, in 
general, the indices affected by prominent sky-lines are poorly fit; 
they were excluded from the fits, but the deviation parameter was 
still computed. For the galaxies with very high central CN indices 
(e.g. central region of N7490), those indices were poorly fit in all 
cases. 

While no single case seems to be ideal for all situations, we 
take the Case Z fits as our best age and metallicity estimates based 
on Lick indices as, in most cases, these appear to be the most stable 
and are the most desirable in terms of excluding highly deviant in- 
dices from the fits. Still, deriving SP parameters from Lick indices 
for star forming galaxies is clearly not a secure method. With our 
radially resolved spectra, we can look at radial trends and identify 
the bins that seem most deviant, but if we should consider only one 
bin per galaxy, we would not be able to make such judgments. Thus, 
let us now consider SP fits to the full spectrum in hopes of finding 
a more stable method for the measurement of age and metallicity 
in star forming galaxies. 



B3 Age and Metallicity Fits from Full Spectra Single SSP 
Fits 

We now examine full-spectrum fits, considering here just the best- 
fit single model SSP. Again we use the BC03 models with all 
model and galaxy spectra convolved to the same effective resolu- 
tion. While fitting the full spectrum makes use of all the information 
in the spectrum, there are some drawbacks as well. Most notably, 
the wider wavelength-baseline reinforces sensitivity to dust extinc- 
tion and relative fluxing errors. However, in this case we can at- 
tempt to model the dust content by incorporating a prescription for 
dust in the models. The fits performed here to single SSP templates 
are done in a very similar manner to the full spectrum fits described 
in §3, and we refer the reader to that section for full details of the 
dust models and the code. Here, the number of templates in each fit 
is M = 1, and the best fit is determined by searching the x 2 -space of 
the fits to each of the 70 individual SSP templates shown in Fig. [2] 
Additionally, for each SSP template, we first compute the x 2 of the 
best-fit dust-free model. We then consider the same SSP, but with 
increasing amounts of dust, as described in £13.4.11 At each dust 
level (increasing by increments of t v = 0.5 at each step), if the \ 2 
decreases from the previous fit, the dust content is increased and 
the SSP refit. This procedure is repeated until we reach a minimum 
in the x 2 -space and only the lowest-^ 2 SSP+dust combination for 
each template is kept for consideration in the final \ 2 grid search. 
We imposed a maximum dust content of fy = 8.5, but this limit was 
never formally imposed, having always reached a minimum x 2 at 
a dust level lower than the limit. 

In Fig. lB3l we show two examples of our SSP fits to the central 
spectra of U2124 (left) and N7495 (right). The central spectrum of 
U2124 shows no evidence for emission when examined by-eye (see 
Fig. lAlt . and the SED is generally well represented by a very old 
(20 Gyr) but low-metallicity (Z = 0.004) SSP with a small amount 
of dust reddening (fy =0.5). The spectrum of N7495, on the other 
hand, shows significant emission, even in the higher-order Balmer 
lines. While these do get masked out by our cr-clipping procedure, 
this spectrum is clearly not well represented by a single SSP. The 
entire spectrum blueward of ~ 4700 A is masked out in the fit, there 
being no SSP that can simultaneously match the low- and high- 
frequency characteristics of the spectrum. 

Further insight into the single SSP fits is revealed by examina- 
tion of the x 2 distribution. In Fig. lB4l we show the corresponding 
X 2 distributions to the fits shown in Fig. IB3I For U2124 (left), the 
lowest x 2 is at the 20 Gyr Z = 0.004 SSP (indicated by the black 
cross). However, a clear minimum in the x 2 distribution has not ac- 
tually been reached. Since 20 Gyr is the maximum age of the SSP 
models, we cannot say whether the x 2 would continue to decrease 
for even older SSP age It is also clear that the x 2 f° r me 2 Gyr 
Z = 0.05 SSP is not all that different from the absolute minimum. 
This is a clear example of the age/metallicity degeneracy in stellar 
populations whereby an old/metal-poor SSP is as good a fit as a 
young/metal-rich SSP. Comparing the results to the full synthesis 
fits (Fig. [3), we see that they provide much better overall matches 
to the SEDs. Looking at the corresponding SFH for r = 0" (Fig.|7]l, 
we see that there is about an equal contribution in terms of I 7 -band 
light-weight from an old/metal-poor SSP, and old/metal-rich SSP, 
and a young/metal-rich SSP, indicating that this spectrum is not 
well-represented by a single SSP. On the other hand, the x 2 distri- 

8 Note that we are not concerned here with model ages that are older than 
that of the Universe. Model ages are not precisely calibrated and, as such, 
we are primarily concerned with relative trends. 
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Figure Bl. Results for age and metallicity fits from Lick indices for U2124. See text for plot description. 



R = 0.0" R = 4.2" R=18.1" 

Age (Gyr) 12.7 0.011.812 5 0.0 20.0 1.420.0 12.519.5 10.0 1.0 9.310.0 7.2 
l°gio( z /Z Q ) -0.7-0.4-0.7-0.8 0.0 - 1,0-0. B- 1.0- 1.2-1.0 -1.0 0.0-1.0-1.2-0.8 



[MgFe]' ; -B 4 



H<5 P 
H7. 

HS, 
TiO, 
TiO, 
NaD 
Fe5782 
Fe5709 
Fe5406 
Fe5335 
Fe5270 
Mgb 

Mg, 
Fe5015 
H/S 
Fe4668 
Fe4531 
Ca4455 
Fe4383 
04300 
Ca4227 
CN, 
CN. 



S> B- - 

b - 

-T_£ -D-p 



^ 



£*A 

-js 

— -i -at 

; *H ~ B — 

Bts- -b|- -a — 
«l-b 

i*. 

f- -U*B 

i» 

-( - as 

^-B 

— ^ — : 

^ 



a-; 



- e-a: 



* - - 



-4- 



• - -HFJ- 

m 

-a -9-1 



-* 

« 

EjXi 

it 



<*B 



-B^ 

m 

jw 

<M 

a 

-w 

-H^L 

B- ft' 

h B-oa- 

rjss-j 

- s — -«S-i 

^ 

J* 

* 

■» 

bX'^ 



: -a. 

; - 



[MgFe]' 
<Fe> 
H7 F 
Hi r 
"7. 
H5, 
TiO, 
TiO, 
NaD 
Fe5782 
Fe5709 
Fe5406 
Fe5335 
Fe5270 
Mgb 

Mg, 

Fe5015 

Fe4668 

Fe453 1 
Ca4455 
Fe4383 
G4300 
Ca4227 
CN 2 
CN, 



-20-10 10 20 -20-10 10 20 -20-10 10 20 

Deviation ( \ = (data^fit,)/*,) N7495 



Rad (") 
s r e 

0.5 I 1 4 5 



0.5 


-0.5 
0.5 

, 
-0.5 
' -1 
-1.5 
20 



-2.4 



* -P~0 — - * 

b. * ,o 5.-.-..#.^©t.-; - x -.. ; Z j : 



s. 



mm®* 



H -h 



J£i_ 



-0.4 



0.4 



! -1.6 -1.2 -0.8 

log 10 (Rad/R d ) 
A: all well-measured indices 

B: all - (Balmer, Fe5015, it CN,) Z: Dev |data,-model,|/«, < 1.0 

C: all- (HP it Fe5015) N74g5 

D: non [a/Fe] affected 



Figure B2. Same as Fig.|BT]except for N7495. 



bution for N7495 (Fig. IB4I right) reached a minimum within the 
SSP parameter space. However, the minimum is not well defined 
and a broad range of age/Z/f combinations are essentially equally 
well fit. Again, from the SFH plot (Fig. [7] right), we see that the 
central spectrum of N7495 has significant contributions from both 
very old and very young SSPs, which cannot be distinguished in 
a single SSP fit. Of course, there are examples of single SSP fits 
where a clear and unambiguous minimum in the age/Z/f \ 2 -space 
exists, but these are the exception for our spectra of star-forming 
spiral galaxies. 



APPENDIX C: COMPARISON OF DIFFERENT FITTING 
TECHNIQUES 

It is already quite clear that the only truly representative fitting tech- 
nique to derive reliable relative SP ages and metallicities for the 
integrated optical spectra of star forming spiral galaxies is the full 
synthesis technique presented in i]3.1l However, it is still useful to 
contrast the results from all 3 techniques to examine any clear rela- 
tive trends, thus allowing the results from each method to be inter- 
preted in the context of the other techniques. In Fig. lCll we compare 
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Figure B3. Examples of single SSP fits to central spectra of U2124 (left) and the emission-line dominant N7495 (right). Colours, lines, shading, panels, and 
labels are as in Fig.|3]butthe ages and metallicities are now SSP-equivalent values. To be compared with full synthesis fits in Fig. [3] 
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Figure B4. Examples of \ 2 distributions for single SSP fits to the central spectra of U2124 (left) and the emission-line dominant N7495 (right). Colours and 
point types indicate metallicity as labeled in the upper right corner of the left panel. Point size is proportional to the preferred fy of the fit to a given SSP. The 
lowest x 2 SSP is indicated by the black cross. 



the derived ages (top panels) and metallicities (bottom panels) be- 
tween the single SSP vs. full synthesis (FS) fits (left panels), the 
single SSP vs. Case Z (Ind "Z") index fits (middle panels), and the 
Ind "Z" vs. FS fits (right panels). All quantities are light-weighted, 
and note that we label the FS parameters as average quantities. Note 
also that the single SSP ages and metallicities are stochastically 
sampled in that they can only be exactly equal to one of the 70 
SSP templates considered, whereas the fits bases on Lick indices 
are interpolated between the model grid space, and the FS fits are 
weighted averages of all SSPs contributing to the fit. 

There is a reasonable correlation between the SSP and FS ages 
(upper left panel) but the slope is far from unity (dotted line). This is 
expected due to the difference between SSP-equivalent ages from 
those derived from the average of a linear combination of SSPs. 
The difference and scatter become more severe the younger the 



age, and the SSP ages reach much younger light-weighted values 
than those of the FS fits. This is once again due to the fact that 
the youngest stars are bright and the single SSP fits are sensitive 
to the stars contributing most dominantly in flux to the integrated 
spectrum, whereas the FS fits are sensitive to any significant con- 
tribution from an underlying old population. On the other hand, 
there is no correlation between the SSP vs. Ind "Z" fits (top middle 
panel). While these two fits both provide SSP-equivalent age esti- 
mates, their primary source of age discrimination is quite different. 
In the full spectrum SSP fits, it comes largely from the overall shape 
of the SED. For the index fits it comes from the age discriminating 
power of the indices that do not get clipped from the fit, which may 
be different for each fit, and are entirely insensitive to the over- 
all SED colour. Finally, the Ind "Z" vs. FS comparison (top right 
panel) shows a weak correlation with large scatter. The spread in 
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Figure CI. Comparison of light-weighted ages and metallicities derived from the three different fitting techniques: single SSP vs. Full Synthesis (FS) [left 
panels], single SSP vs. Case Z Lick index (Ind "Z") [middle panels], and Ind "Z" vs. FS [right panels]. Point types and colours are coded to a given galaxy, as 
labeled in the upper left panel. The dotted line in each panel marks the one-to-one relation. 
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Figure C2. Comparison of the radial gradients from the three different fitting techniques. Panels and symbols are as in Fig. [8] Different line types mark results 
from the 3 fitting techniques, as indicated in the top right of each panel. 



the age predictions is similar between these two methods, indicat- 
ing that neither is overly sensitive to just the most recent episode of 
SF, but the slope of the correlation is not unity, with the FS fits tend- 
ing to older ages. As expected, in all comparisons, the best match is 
found for the oldest ages, corresponding to those spectra for which 
there is no evidence of emission lines (e.g. U2124 and N7490). 

The situation for metallicity is quite similar, with a weak and 
less than unity correlation between the SSP and FS fits (bottom 



left), essentially no correlation for the SSP and Ind "Z" fits (bottom 
middle), and a stronger non-linear relation for the Ind "Z" and FS 
fits (bottom right), with the FS fits tending to higher metallicities. 

Given the weak to null correlations between the derived SP 
parameters between the three methods, it is also important to in- 
vestigate how the derived gradients would be affected. In Fig. IC2I 
we show the radial variations in age and Z for the three different 
fits. The dashed-dotted lines are for the Ind "Z" fits and have er- 
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ratic profiles that can jump drastically from one radial bin to the 
next. These jumps often land at the SP model extremes, indicating 
that they are poorly determined. The metallicity profiles seem to be 
much more stable that the age profiles, and as we saw in Fig. IC1I 
they are biased to much lower values than the FS fits (solid lines). 
The SSP profiles (dashed lines), on the other hand, are much more 
stable from one radial bin to the next. For the relatively emission- 
line free U2124 (left), the age profile has a similar overall form 
to that of the FS fits, but because they are insensitive to SPs pro- 
viding small relative contributions to the SED flux, the SSP bulge 
looks much older than is indicated by the average age from the 
FS fits. Moreover, the metallicity profile is constant at a sub-solar 
Z = 0.004, whereas both other fits indicate radial variation in Z, 
with higher metallicity towards the center. This is another manifes- 
tation of the age/Z degeneracy, particularly in the context of SEDs 
that are made up of SPs of a range of ages and Z's and are thus 
poorly represented by a single SSP. 

The emission-line dominated N7495 (right) shows a slightly 
different set of pathologies. The Ind "Z" age profile is again erratic, 
jumping from very old to very young values at consecutive radial 
bins. The single SSP fit profile is flat, but saturated at the lowest 
SSP age, being only sensitive to the significant amount of current 
SF and blind to any underlying population. The Ind "Z" metallicity 
profile again shows a similar trend as the FS fits, but at systemat- 
ically lower-Z. The single SSP fit profile here is erratic, jumping 
between high and low-Z extremes. 

Finally, the radial profiles in the amount of dust extinction in- 
ferred form the FS and SSP methods is shown in the bottom panels. 
Both are in agreement for U2124 which is consistent with very lit- 
tle or no dust at all radii. The profile shapes for N7495 are similar, 
but the SSP derived dust levels are much higher. The higher values 
of f in the single SSP fits could in part be absorbing the redness of 
the underlying old SP that the single young SSP does not accom- 
modate. 

To summarize, there is no clear correlation between the three 
fitting techniques for any spectrum with significant contributions 
to the light from multiple SSPs. The SSP-equivalent parameters 
derived from either the index-based or single SSP fits reveal erratic 
behaviors in the radial profiles, which are very sensitive to different 
sets of characteristics of the spectrum. Only the individual fits and 
radial profiles from the full spectrum population synthesis fits yield 
representative and stable results and are suitable for assessing rela- 
tive trends in age and metallicity. Therefore, the analysis and results 
presented in the main paper are based solely on the full spectrum 
synthesis fits. 
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